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PLASMA JET PIERCING OF MAGNETIC FIELDS AND ENTROPY TRAPPING 
INTO A CONSERVATIVE SYSTEM* 


James L. Tuck 
Los Alamos Scientific Laboratory, Los Alamos, New Mexico 
(Received August 17, 1959) 


To obtain a thermonuclear plasma inside a 
magnetic confinement system, we have either 
(a) to create the plasma in situ from cold gas by 
ionization and heating or (b) build up the plasma 
piecemeal by injection from some accelerator or 
hydromagnetic gun. 

Both of these methods have characteristic re- 
quirements: For (a) as in the pinch effect, great 
input power is needed in order to get through the 
strongly radiating partly ionized plasma phase 
and after that it is difficult to raise temperatures 
of the ions much above a kilovolt by Joule heat- 
ing on account of the high electrical conductivity 
of the plasma and the slow exchange rate between 
the electrons which receive the Joule heating and 
the ions which it is desired to heat. For (b) 
single particles cannot of course, be trapped 
into a static magnetic field and some method of 
circumventing the conservative property has to 
be found. Trapping by means of a rapidly opera- 
ting magnetic trap door has been discussed by 
Post. Another nonadiabatic process can be the 
breakup of an energetic molecular ion such as 
D,*-D* + Dt +e by collision with an internal arc 
column as in the DCX experiment (Luce) or by 
collision with residual gas as in the OGRA ex- 
periment (Golovin). Energetic neutral D, can in 
principle be used instead of D,*. Still another 


process can be by interparticle collisions, as 

for example the injection of plasmoids by Bostick. 
Such plasmoids soak into the magnetic field and 
after developing the appropriate Hall polarization 
can continue across it. In this note, we propose 

4 new method of achieving the nonadiabatic pro- 





cess which may prove very advantageous. 

Consider a collisionless plasma of pressure 
p, confined on all sides by a magnetic field of 
magnitude (87p)”?. Let the plasma now assume 
the improbable but possible configuration in 
which the particles all are directed at some 
small region on the wall. We see that the parti- 
cle pressure there becomes greater than the 
magnetic field pressure, so that the field can be 
forced aside letting plasma escape. If now we 
review these events in reverse time sequence, 

a method emerges for the trapping of plasma 
into a magnetic field. 

Consider the interaction of a broad parallel jet 
of cold plasma, density p,, velocity v, in the x 
direction, electron mass m, ion mass M, inci- 
dent on a magnetic field of magnitude B of large 
radial extent. The velocity v, is large, such that 
the cross section for ion-electron collisions with 
relative velocity v, is negligible (no collisions). 
These conditions resemble those for the plasma- 
field sheath described in the infinite-conductivity 
theory of the pinch (Rosenbluth and Garwin)' ex- 
cept that the particles here may be incident at 
angles other than normal to the magnetic field 
direction. According to this model, charge sep- 
aration occurs, setting up an electric field which 
retards the ions and speeds up the electrons. 
The deflection of the ions is predominantly elec- 
trostatic, so that they are brought nearly to rest; 
the electrons reach an energy close to the orig- 
inal ion energy in the jet, and provide most of 
the electric current which screens off the mag- 
netic field. The thickness of the sheath over 
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which the magnetic field and plasma density drop 
reciprocally to zero is given by d=v,c(mM)”?/eB. 
Particles enter and leave the sheath at the same 
angle. 

We choose the stagnation pressure of the jet, 
PoVo", to be greater than the magnetic field pres- 
sure so that yielding of the field occurs. (For 
simplicity in the order-of- magnitude discussion 
which follows, the field is chosen of large extent 
so that the magnetic field cannot pile up to larger 
values.) Let us consider the situation where the 
jet has pierced the magnetic field for some dis- 
tance. We shall assume that the process of re- 
flection off the curved end of the cavity can be 
reasonably approximated by uniform (cosine law) 
scattering from a plane end. 

Let the rate of elongation of the cavity (“tip 
speed”) be v. Writing the pressure balance equa- 
tion for the end [remembering that for cosine 
scattering the mean backward particle velocity 
is one-half the particle speed (v, - v) relative to 
the wall], then 


$ po(v, - v)? = B?/8n, 


or 


tip speed i ; 
~ jet speed = Poo?} — b- 8/85 


B? = 12npvg? = 7.7107 vE, (1) 


where j =deuteron equivalent current density in 
amp/cm? and E = deuteron energy in ev. Similar 
equations were derived for the neutral stream 
magnetic storm theory of Chapman and Ferraro’ 
in 1930. We see that a critical value B, of the 
magnetic field will bring « and the penetration 
velocity to zero. Larger B give a negative e, 
corresponding to extrusion of a previously pene- 
trated beam. Using mass conservation and let- 
ting p be the internal density after scattering, 
we next write 


PoUo= ple (vy - v) - v]+pv, 
or 
Pp =2pov,/(v, - v). (2) 
The radial pressure on the walls, holding apart 
the magnetic field, using this density and the 
mean radial velocity $(v,- v), turns out to be 
PoVo(¥_ - v) which is only 2v,/3(v, - v) of the re- 


quired value of B?/8z. To meet this difficulty, 
(a) the cavity could expand leaving a neck at 
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entry —this neck reflects back some of the out- 
wardly escaping particles raising p to the value 
required to keep the walls from collapsing, or 
(b) the walls could contract to the point where 
they intercept some of the incoming beam, again 
raising the wall pressure. Resorting to the anal- 
ogy of the cavitation of water by a jet of air, we 
find both cases reported in the literature of 
hydrodynamics together with a third possibility- 
Rayleigh-Taylor instability with mixing. We may 
expect the latter to occur in magnetic field pierc- 
ing especially when the magnetic field direction 
is transverse to the jet direction, since in such 
a case the magnetic field-plasma boundary is 
necessarily curved in the unstable sense.’ Such 
instabilities will incidentally trap plasma and 
may be responsible for the trapping reported by 
Coensgen and Ford.‘ For injection in the axial 
direction, the possibility of hydromagnetic sta- 
bility appears to be greater; Helmholtz insta- 
bility is minimized by the highly supersonic 





state of the jet.° It might be argued that one kind 
of trapping is as good as another, and thatthis | 
process can be welcomed; as will appear how-__ | 
ever, we require the plasma to attain the maxi- 
mum possible density with a minimum of inter- 
mixed magnetic field. 

In order to exploit the trapping possibilities, 
consider a beam entering a magnetic field such 
as a picket fence or cusped geometry® which 
after the initial increase to some value less than 
B,, falls again. The beam can pierce the barrier 
forming a neck but in doing so must raise its 
pressure to a value larger than the field pres- 
sure further on. Consequently the jet cavity can 
expand against the field further down stream, 
forming a balloon (Fig. 1). Any jet which makes 
a clear traverse of the peak field without deflec- 
tion, can be scattered off a B>B, “anvil” field 
further on, though this is not essential to the 
argument since scattering at the entry reduces 
the beam piercing power. The size of the plasma 
balloon will be governed by the balance between 
input flow and loss by back flow and diffusion 
into the magnetic field. 

If the jet is now stopped, pressures fall and as 
the magnetic walls close in, the neck closes in 
faster, pinching off the balloon and leaving a 
plasma bubble trapped inside. The picket fence 
seems particularly well adapted to this method, 
since the field has the appropriate valley in the 
center of the machine, though the method seems 
adaptable with trivial modifications to Stellarator 
and mirror machine geometries.’ Two types of 
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out- for which no further heating is required, and we 
value see that for E=90 kev, J~100 amperes. This is 
_ or a current outside present accelerator practice 
re B<B B>B by a factor of 100 but should not be too difficult, 
again Cc c if one uses a condenser power supply. More 

> anal- immediately available seems to be a hydromag- 
r, we netically accelerated plasma jet. Taking tenta- 
f tive estimates of particle density 10** cm“, 

lity - — V,= 10" cm/sec for quantities already achieved, ® 
/e may 2A we find the critical magnetic field for piercing 
Pierc- to be a sizable 6500 gauss. This hydromagnetic 
ction gun already produces a volume of jet adequate 
such for inflating several liters to a density of 10°° 

is cm~*. A 30-fold increase of the present velocity 
Such would bring the plasma to the optimum tempera- 
ind ture for a DT thermonuclear reactor. 

ed by It is a pleasure to acknowledge helpful scienti- 
xial FIG. 1. Jet piercing into a picket fence magnetic fic discussions with C. L. Longmire, J. Marshall, 
sta- field . F. L. Ribe, and M. Rosenbluth. 

ita - 

Cc operation come to mind: (1) pulsed inflation . 

e kind | whereby the plasma bubble pinched off during ‘Work performed under the auspices of the U.S. 





this | the off-beam period by the collapsing walls of 


ow- | the entry channel is left to react inside the ma- 

naxi- | chine; (2) steady operation whereby a large 

iter- volume is kept full against losses by mixing and 
back-streaming. In the latter case, the prime 

ies, consideration for a thermonuclear reactor be- 

such comes the chance of an ion reacting compared 

h with its chance of finding the exit channel and 


sthan escaping. At first sight it might appear that this 
yarrier ratio might be made arbitrarily advantageous by 


its reducing the jet to needle-like dimensions as 
‘eS- from an accelerator. However, there is a limit 
ry can | tothe narrowness of the beam set by (1) break- 
m, down of the infinite plane sheath approximation 


nakes | and (2) large pressure balance errors when the 
eflec- | sheath thickness d becomes comparable with the 
field jet diameter; and furthermore boundary layer 

he instability may erode the jet during the penetra- 
uces tion, limiting the ratio of the penetration depth 
plasma | of the unscattered beam to the jet diameter to 
ween | SOme constant value as in the hydrodynamic ana- 
on log. 

The expression for the total current carried by 
and as } one kind of particle in a deuterium jet of cross- 
2s in } Sectional area ad’, using the current density 

a from Eq. (1), turns out to be (independently of B) 
wee 1=1.35x10"5q@*mc72VE, (E in ev). 

nthe } We do not yet know how low a canbe. Making a 
seems | Teasonable choice of a = 10, we have 

larator 


wr 1=0.3VE amperes. 


We are naturally interested in injecting plasmas 
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"The picket fence geometry, long known to have great 
hydromagnetic stability, has received little experi- 
mental study mainly because of its obvious disadvantage 
with respect to the other main geometries, pinch, 
Stellarator, mirror machine, of leakier confinement. 
The prediction of strong cyclotron radiation from a hot 
plasma containing a magnetic field by B. A. Trubnikov 
and V. S. Kudryavtsev [Proceedings of the Second Uni- 
ted Nations International Conference on the Peaceful 
Uses of Atomic Energy, Geneva, 1958 (United Nations, 
Geneva, 1959), Vol. 31, p.93], may make it of pri- 
mary importance to avoid this radiation by using plas- 
mas containing zero or little magnetic field. Only the 
picket fence has enough hydromagnetic stability to con- 
fine such a plasma. The other devices are theoretically 
limited to 8 (= plasma pressure/magnetic field pressure) 
less than about 0.1. 

8J. Marshall (to be published). 
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FERMI ACCELERATION OF ELECTRONS IN THE OUTER VAN ALLEN BELT* e 


J. A. Crawford 
Berkeley Astronomical Department, University of California, Berxeley, California 
(Received September 8, 1959) 


The persistence of the inner Van Allen belt is 
well understood to be due to the constancy of the 
strong geomagnetic field close to the earth. The 
weak field in the outer belt, however, is violently 
disturbed about once a month by a large solar 
storm cloud ejected by the sun in association with 
a flare. One can show that this solar wind should 
in one way or another wipe out any radiation 
originally present in the outer belt. In order 
therefore to account for the presence of this in- 
tense radiation it has been suggested that it is 
originally trapped in the cloud and deposited in 
the earth’s field during the geomagnetic storm.’ 
It will be shown here that this assumption is 
unnecessary since the radiation in the outer belt 
may be accounted for by a local Fermi accelera- 
tion created by the interaction between the storm 
and the geomagnetic field. Parker has shown 
that this interaction gives rise to instabilities 
which are most pronounced at the smaller wave- 
lengths. Thus, in the region where the energy 
density of the storm and the earth’s field are 
equal, turbulence should consist mostly of eddies 
of the smallest possible wavelength, i.e., a pro- 
ton Larmor circumference.’ While the storm is 
under way, we may expect that closed magnetic 
loops will frequently be broken off by Sweet’s 
mechanism.’ The electrons trapped in these 
loops will be able to undergo a Fermi accelera- 
tion without being absorbed in the earth’s atmos- 
phere as would quickly happen if they were trapped 
on lines connected to the earth. It is this circum- 
stance which permits them to reach very high 
energies. Protons, on the other hand, are not 
effectively trapped in the loops because of their 
twisty structure: a twist should have a length of 
the order of a proton Larmor radius and there- 
fore scatter protons magnetically. Thus we do 
not expect the proton energies to reach the Mev 
range through a local Fermi acceleration. In- 
stead, protons will diffuse and be gradually ex- 
changed with others, lobbed upward from the 
top of the atmosphere.* The loops remaining 
after the storm has passed will be enmeshed in 
the geomagnetic field and will eventually decay, 
depositing their electron radiation into this field 
where it will remain trapped by the mirror mech- 
anism until the next storm destroys it by dump- 
ing it upon the earth through a Fermi accelera- 
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tion along the geomagnetic lines of force. 

In order to show that this mechanism can real}; 
work, one must (1) estimate the magnitude of the 
Fermi acceleration and (2) take proper account 
of the collisions between electrons which tend to 
work against it. This has been done by setting up 
and solving a suitable Fokker-Planck equation, 

(1) For a proper treatment of Fermi accelera- ) 
tion one is required in general to take account of 
the reaction stress of the flux of accelerated 
particles upon the plasma adhering to the fluctua- 
ting lines of force.* This makes it possible to 
obtain an equilibrium energy spectrum for the 
fast particles. In the present instance we may 
assume that this equilibrium is frozen in when 
the Fermi acceleration produced by the departing 
cloud has become so weak that the acceleration 
time scale for typical Van Allen electrons is of 
the order of the decay time of the storm. We t 
may take this as about a half-day (4x10* sec). ? ‘ 
In this time a slow electron reaches on the aver- tl 

tl 





— a. ae. oe ee 


age a speed of 10'° cm/sec (~30 kev), correspon 
ing to a Fermi acceleration, ¢ =2.5 10° cm/sec’ 

(2) The equation of motion for an electron of 
momentum p in a steeply helical path around a 
moving line of force is h 


. ({B =B ' 
Q =a-| —-V— 

3(5-#), 
@ being the electric drift velocity and B the mag- 
netic field. Defining E =p?/2m, which is the 


energy at moderate velocities, we find, if the 
magnetic line moves in a random manner, 


f, = lim Sa82 -26(22 ) ; 


At~0 At v : 


p =Qp, 


——— 
~~ 


/ 2 2 
6 = Me ( (AE )*) mx» § 
‘es At~0 At v 





where v is the electron velocity and x = -(2) 


represents the reaction of the fast electron flux z 
on the plasma. @¢ is given by z 
oc \ 
¢ -4f dz(Q(s, t)2(s +z, t)). 
—_ G 
Here s is the electron coordinate along the mag- fh 
u 


netic line at the instant ¢. At nonrelativistic 
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energies, 





lim (Av) 


at—o 4 — 


the mean acceleration. The parameters f,, f, 
determine the Fokker-Planck equation for the 
distribution of E. However, we have neglected 
collisions between electrons. These may be 
taken into account by adding to f, and f, the cor- 
respondingly defined parameters arising from 
collisions. These are well known® and, in the 

£ range where both contributions are compara- 
ble, we may write 


f, =2E[2(¢/v) -x]-Z/E™, 
fo =8E* > /v +2ZkT /E™?, 
Z =4ne‘n InA/(2m)”, (1) 


where A is the Debye shielding parameter. At 
higher energies the effect of collisions may be 
neglected in comparison with Fermi accelera- 
tion; on the other hand, at smaller energies the 
values of f,, f. given by (1) are incorrect; never- 
theless they may be used in computing the steady- 
state solution of the Fokker-Planck equation since 
they give the correct (Boltzmann) distribution in 
this range. The Fokker-Planck equation, 

OF(E,t) 3 


ctl th ate 
at ap? (2heF) op (fF): 


has the steady solution 
F(E) « f, ‘exp f(2f,/f.)dE. 


The integral can be evaluated in overlapping 
sections, yielding the approximately normalized 
result 





2 (B\” [1+66"})2 tan“!y£ 
FI) =-p;—(— -— 
we; eran) 1+ (EP W/me * [ ykT 


Ww »\77xXe/@ 
| 3 (: +=) ; 
mc c 


y =(2¢/ZRT)”*(2m)™. (2) 





The nonrelativistic approximation to this spec- 
trum, in the range y""«<E <m_e’, yields the 
exponential velocity distribution 


~(1/2ykT) ~xv/> 


G(v) = 2(2m /nkT)“?(ykT) “e (3) 


This may be used to estimate the mean energy 
flux density per electron, since the contribution 


of relativistic electrons is not a major one: 
J 


_—~ 


dv vmv’G(v) 


yates) 
y\r kT x 


At not too great geocentric distances, the elec- 
tron density » of the hydrogen plasma left behind 
by the departing storm may be obtained by writing 


nMU? /2 ~ B? /8r, (4) 


where M is the proton mass. If the turbulent 
velocity is equated to the storm velocity, taken 
to be 500 km/sec, then, in the center of the outer 
belt where B~10~* gauss, we obtain n =2 x10° 
electrons/cm*. The effective electron tempera- 
ture left behind by the storm is obtained by 
writing 

3kT = $mU?, 


yielding kT ~0.5 ev. Thus Z can be evaluated. 

@ was previously estimated to be 2.5 10° cm/ 
sec’. Substituting this value in the expression 
for y, we find y"'~10 ev. This value is of course 
not reliable within at least a factor of 2. Within 
this range the observed value of J (2x10" ev/cm? 
sec sterad)’ can readily be obtained for any rea- 
sonable value (between 1 and 10) of the numeri- 
cal parameter xc/®. 

The presence of a gap between the two zones 
appears to be well accounted for by the sweeping 
action of the storm, regardless of whether the 
outer Van Allen radiation is brought from the sun 
or locally accelerated: the region just below the 
level to which the loops penetrate will be depleted 
of radiation by a Fermi acceleration, occurring 
along the geomagnetic lines, which dumps the 
trapped particles upon the earth. This region 
may be identified with the observed gap. The de- 
crease in intensity at large geocentric distances, 
on the other hand, may be due to a great increase 
in the reaction parameter x, as one goes farther 
out. In such regions, where the density is low 
and collisions relatively unimportant, all the 
electrons would like to be accelerated to high 
energies. Since the energy density available is 
much too small for this to happen, it appears 
that a large reaction parameter must be built up 
to prevent it. This will make the high-energy 
spectrum very steep. Thus it seems that, in 
order to produce electron radiation locally in the 
outer Van Allen belt, a Fermi acceleration is 
needed which is neither too weak nor too strong. 
The breakdown of the longitudinal invariant of the 
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electron’s motion around a loop may also fail at 
large geocentric distances, with a consequent 
sharp reduction of the Fermi acceleration and 
therefore of the intensity. 

The writer wishes to thank Professor K. M. 
Watson, Professor L. G. Henyey, Dr. S. Trehan, 
and Dr. W. Hess for their encouragement and for 
some very stimulating discussions. 
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ION TEMPERATURE IN SCYLLA, AS DETERMINED FROM THE REACTION D(d, p)T* 


D. E. Nagle, W. E. Quinn, W. B. Riesenfeld, and W. Leland 
Los Alamos Scientific Laboratory, University of California, Los Alamos, New Mexico 
(Received September 14, 1959) 


Scylla, the Los Alamos experiment on fast mag- 
netic compression of a plasma, has been shown 
to produce about 10’ neutrons per pulse from the 
reaction D(d,n)He*.’ An equivalent ion tempera- 
ture of about 1 kev corresponding to this yield 
can be computed, using our measured values of 
the plasma density (6x 10"*/cm*), volume of the 
emitting region (3 cm*), and mean duration of the 
neutron pulse (0.9 usec). As has often been re- 
marked, however, temperatures based on the 
reaction rate alone deserve little credence until 
the rate is so high as to be inexplicable on any 
other assumption than the thermal one. The pre- 
sent experiment was designed to measure ina 
more direct way the ion temperature. By study- 
ing the velocity spectrum of products of the reac- 
tion D(d,p)T we measure the mean square veloc- 
ity of the center of mass of the d-d system, and 
hence the temperature of the plasma. We find in 
this way T=1.3 kev. The spectrum obtained, 
coupled with the result that the neutrons emitted 
in the radial direction showed no Doppler shift 
on reversing the applied field,’ argues strongly 
against many of the mechanisms for spurious 
neutron production which have been proposed. 
Radial hydromagnetic shocks having speeds of 15 
cm/sec are observed at the time of zero ex- 
ternal magnetic field. It is easy to account for 
the measured final ion temperature with a model 
in which the plasma is adiabatically compressed 
after being preheated by the shocks. The assump- 
tion that the plasma is completely diamagnetic 
leads to the numerical result that the particle 
pressure is in balance with the confining magnetic 
field. 
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For a plasma of temperature T, the specific 
reaction rate is (ov) Av’ namely the reaction 
cross section times the relative velocity aver- 
aged over a Maxwell distribution. The velocity 
spectrum of the protons depends on v,, the mini- 


mum velocity of emission of a proton in the cen- ) 


ter-of-mass system, on kT, on the particle 
masses Mp, Ma, and ™ 45 and on the Q value of 
the reaction. We derive for the number of pro- 
tons per unit energy interval 


m v2 
dN/dE = (4) (n ov), /™ 5%o) 


4nkT 
m 
xexp} if (v - | 


oo wasps: 
="ol"* 30° dinkT 


For protons 
2 = 
m 3% /2 =m ,Q/(m,+m ), 


and for tritons the subscripts p and ¢ are inter- 
changed. By measuring the width of the velocity 
distribution of the products we obtain a value for 
ae 

The axial direction was chosen because no Wil- 
dows in the discharge tube are required, and 
because the interpretation is simpler. The rest 
of the geometry was fixed by the parameters of 
a magnetic spectrograph which was available. 
Figure 1 shows in cross section the experimental 
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FIG. 1. Experimental arrangement. The discharge 
occurs within the alumina tube at the center of the 
coil. The dashed lines represent trajectories of the 
protons and the tritons from the d-d reaction. 


arrangement. The Scylla discharge takes place 
inside the alumina tube. The spectrograph en- 
trance slit is placed as close as possible to the 
discharge, namely 10 cm. The slit, if too near, 
is damaged by heat and wall material which 

strike it during late stages of the discharge. The 
slit forms the object for the ion-optical system, 
and the dispersed image is formed on the nuclear 
emulsion. 

Two runs were made with the spectrograph and 
charged particles from Scylla. A thin cover was 
used to protect the emulsion from scattered light. 
In the first exposure the cover was1.0 mg/cm? 
of (aluminized) polyester film; this effectively 
stopped the tritons, but left the protons with 
about 67 microns residual range. In the second 
exposure the cover was 0.55 mg/cm? of nickel; 
this allowed the tritons to enter the emulsion with 
aresidual range of about 5 y. The first exposure 
was 1617 discharges of Scylla, the second 1888, 
and the corresponding total numbers of neutrons 
emitted, as measured by a silver counter, were 
1.0x10"° and 1.1x10"°. The plates after develop- 
ment showed tracks which were easily identified 
as the d-d protons or tritons, respectively. The 
density of tracks when plotted vs position re- 
vealed smooth distributions whose maxima came 
at the predicted locations. The width of the 
curve in the least favorable case was three times 
the resolution of the instrument and in the most 
favorable case eighteen times the resolution. 

The resolution and energy calibration were meas- 
ured by means of an exposure to a thin source of 
Pu“ alpha particles. The track location fixes 

the particle velocity, and the number of tracks 

Per unit length along the plate is converted into 
lumber of particles per energy interval coming 
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FIG. 2. Velocity spectra of protons and tritons from 
the d-d reaction in Scylla. The solid and the dashed 
curves are computed with the expressions given in the 
text. The width of the triton curve is eighteen times 
the measured resolution of the spectrograph. 


from reactions in the plasma. The conversion 
takes into account the variations along the plate 
of angle of entry into the emulsion and of lateral 
magnification. Figure 2 shows the reduced dis- 
tribution dN/dE for the second exposure. The 
curves are computed ones using the equation for 
dN/dE for T=1.0 and T=1.5 kev. A temperature 
of T=1.3 kev is in agreement with these data and 
with similar data from the other exposure. The 
agreement of » with predicted values is as fol- 
lows: for the protons 2.410x10° cm/sec pre- 
dicted, observed 2.412 x10° cm/sec, and for the 
tritons 0.8047 x 10° cm/sec predicted, observed 
0.8041x10° cm/sec. This agreement indicates 
that the plasma has negligible mass motions in 
the axial direction. The principal uncertainty in 
the interpretation of the distributions is intro- 
duced by the low-energy tail to the proton distri- 
bution; however, this involves only about 1.5% 
of the reactions. 

A more complete discussion of this experiment 
is being prepared for publication. 

We acknowledge with pleasure the contributions 
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TWO MASER EXPERIMENTS TO TEST GENERAL RELATIVITY 


Huseyin Yilmaz* 
Institute for Advanced Study, Princeton, New Jersey 
(Received August 5, 1959) 


Despite the fact that special relativity is firmly 
established on experimental grounds, the same 
cannot be claimed for the general theory. First 
of all, there are only three experimental tests 
of general relativity and not all of them are 
sufficiently conclusive in favor of the theory. 

(i) The bending of light rays, though correct as 
an order of magnitude, is just outside the experi- 
mental error. (ii) The observed shift of spectral 
lines towards red is qualitatively in agreement 
with the theory but a quantitative agreement is 
still a much discussed and investigated question. 
(iii) In the case of the advance of the perihelion 
of Mercury, the situation is different. Here, we 
have a satisfactory quantitative confirmation. 
Until recently, it was thought that the advance 

of the perihelion of Mars was not correctly pre- 
dicted by the theory.’ In the last few years, 

more elaborate calculations have shown that the 
theory predicts the effect within the observational 
error.” The calculations for the Earth’s perihe- 
lion are at present in progress, but preliminary 
calculations are in agreement with the theory. It 
must be remembered that the perihelion motion 
is a second order effect and its quantitative agree- 
ment lends great support to the theory. Never- 
theless, one cannot say that these facts provide 
sufficient tests for such a fundamental physical 
theory as general relativity. There is definite 
need for further practicable experiments in the 
field of gravitation. This need is being felt all 
the more strongly as time goes on because of 
man’s increasing interest in space and gravita- 
tion. 
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In this note, two closely related but essentially 
different experimental tests of the theory are 
discussed. Only a year ago, it would have been 
impossible to carry them out. Thanks to the 
recent advances in maser techniques, these ex- 
periments can now be done without too much 
difficulty.® 

The first experiment* provides a direct test 
for the principle of equivalence. As is well known, 
according to the principle of equivalence, the ef- 
fects of an external gravitational field in a local 
coordinate system (say a small material box) can 
be eliminated completely by letting the coordinate 
system be free in that gravitational field (freely 
falling box). In such a coordinate system, the 
external gravitational field does not have any 
influence either on the motions of particles or 
on any other physical process whatever. Thus, 
according to the principle of equivalence, such 
a system must locally be equivalent to a Lorentz 
frame and in it, light rays must travel with the 
same velocity in all directions. Now, in the 
gravitational field of the Sun, our Earth may be 
considered as a small box and the above state- 
ment concerning the local propagation of light 
may be tested by comparing the velocity of light 
in the direction of the line joining the Earth and 
the Sun and in a direction perpendicular to it and 
the Earth’s radius (Fig. 1). Nowadays, maser 
techniques are approaching an accuracy 
5c/c =10~ (c is the velocity of light) in com- 
paring the two velocities.* If a discrepancy is 
found to this accuracy, it would imply a devia- 
tion from the principle of equivalence in first 
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FIG. 1. SV is vertical; SM, and SM, are horizontal. 
M, and M, are mirrors and S is the light source. 


order of the dimensionless quantity GM/c?R 
~9x10-°, where M is the mass of the Sun, G the 
gravitational constant, and R is the distance of 
the Earth from the Sun. This dimensionless 
quantity is the expansion parameter involving 

the Sun’s gravitational field. A higher order 
deviation would be presently undetectable. 

The second experiment is similar in appear- 
ance but essentially different in nature. To see 
its function clearly, let us assume that the first 
experiment gave a negative result so that the 
principle of equivalence is confirmed. We can 
then ask the following question: Is the propaga- 
tion of light the same in all directions in a co- 
ordinate system fixed in a gravitational field? 

A laboratory room on our Earth is such a fixed 
coordinate system in the gravitational field of 
the Earth. This question can be answered by 
comparing the velocity of light in the vertical 
and horizontal directions on the Earth (Fig. 2). 
If there is a first order influence due to the 
Earth’s gravitational field, it would be of the 


Yu “2 








FIG. 2. SM, is horizontal and SM, is vertical. M, 
and M, are mirrors and S is the light source. 


order of the dimensionless quantity Gm/c?r 
=7x1071°, where m is the mass of the Earth and 
y its radius. This experiment does not test the 
principle of equivalence. It tests an implicit 
assumption, namely, the local isotropy of the 
space-time continuum. This assumption may 
conveniently be expressed as follows: The grav- 
itational field influences all lengths and length- 
measuring rods in the same way. Also, all 
periodic phenomena and time measuring devices 
are influenced in the same proportion. As a re- 
sult, one never is able to detect locally the effects 
of a gravitational field by kinematical means. 

The second experiment is considerably more 
difficult then the first one because when one 
rotates the maser into the vertical directions, 
stresses and strains are introduced due to the 
Earth’s gravitational field. In the first experi- 
ment, such complications can be avoided by 
doing the experiment during the sunrise or sun- 
set. In this case, rotations are only around the 
vertical axis.°® 

In conclusion, we would like to emphasize that 
the experiments discussed here are essentially 
independent of any theory although it would have 
been difficult to present them without the aid of 
the general theory of relativity of Einstein. 

The author would like to acknowledge a dis- 
cussion with Professor C. H. Townes who made 
many essential suggestions. 
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(1947). Also, private communication with the Naval 
Observatory, Washington, D. C. 

3C. H. Townes and A. L. Schawlow, Phys. Rev. 112, 
1940 (1959). See also J. Cederholm et al., Phys. Rev. 
Letters 1, 342 (1958). 

‘Experiment No. 1 is being considered by C. H. 
Townes and J. Cederholm at the IBM Watson Labora- 
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5Suggested by Professor C. H. Townes. 
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KNIGHT SHIFT IN SUPERCONDUCTORS* 


Paul C. Martin? and Leo P. Kadanofft 
Lyman Laboratory of Physics, Harvard University, Cambridge, Massachusetts 
(Received September 2, 1959) 


Although the theory proposed by Bardeen, 
Cooper, and Schrieffer’ has been strikingly suc- 
cessful in explaining most properties of super- 
conductors, the measurements of electron spin 
susceptibilities are supposed to present it with a 
serious dilemma. The BCS theory introduces 
strong coupling between pairs of particles with 
opposite spin and momentum, and therefore tends 
to prevent the electron spin reorientation neces- 
sary for low-temperature paramagnetism. How- 
ever, observations of the electron spin paramag- 
netism, by the nuclear hyperfine interaction” * 
(the Knight shift) in superconducting samples 
small enough to permit magnetic field penetra- 
tion, indicate a substantial spin susceptibility. At 
absolute zero, the spin susceptibility observed in 
tin and mercury is approximately three quarters 
of the free electron value. 

The purpose of this Letter is to suggest a mech- 
anism for illuminating the Knight shift within the 
BCS framework. This mechanism occurs in con- 
sequence of the relative magnitude at low tem- 
peratures of the typical magnetic field wavelength, 
\ =1/k (approximately 10-5 to 10-® cm), and the 
typical coherence distance, £, (approximately 
10-* to 10° cm). As these numbers indicate, the 
magnetic field contains Fourier components with 
wavelengths small compared to the size of the 
bound pairs. These components tend to polarize 
the pairs and thereby eliminate their compensa- 
tion of the free electron spin susceptibility. In- 
deed, as the typical magnetic field wavevector 
increases, the paramagnetic spin susceptibility 
approaches the Pauli free electron value, y°=3ny?/ 
mv*. (In this expression yu and m are the elec- 
tron moment and mass, » is the electron density, 
and v the Fermi velocity.) 

In the first paper of a series on many-particle 
systems,‘ methods for determining thermal and 
electromagnetic properties of quantum mechani- 
cal systems have been developed. In a second, 
in preparation, these techniques will be shown to 
yield, for superconductors, the BCS solution with 
minor modifications which do not destroy the gap5 
Using this approach, we have investigated the 
electromagnetic properties of superconductors, 
obtaining results which agree with the infrared 
absorption experiments on thin films,® and with 
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the calculations of others.” We have also applied 
these methods to the electron spin paramagnetism 
of a superconductor in a static but spatially var- 
ying magnetic field, H. In sucha field, the spin 
magnetization, M, is characterized by a non- 
local susceptibility, x, which satisfies 


M(r) = Sipe - ¥’)H(r’)dr’. (1) 


The Fourier transform of this spin susceptibility, 
x(k), has been evaluated in terms of the attenua- 
tion parameter, a= 3hvk/A(g), the half energy 
gap, A(g), and the inverse of the absolute tem- 
perature in units of the energy, 8. The result, 
derived by integrating the time-dependent corre- 
lation functions of the BCS model, may be ex- 
pressed as 





1 
R(a, p) =1- 4a tanh(3A) 
1 fa tanh[3gA(14x2)"7] 
pa X, |a+x pBA(1+x 
* 3a x - a-x (14x?) ¥2 » (2) 
0 








where R is the ratio, y(k)/x°. 


In the large-attenuation or Pippard limit (a>>1), 
and at low temperatures, this ratio becomes 


R(a>>1, B—0) =1 - 1?/4a =1 - mr/2E,. (3) 


The experimental results are fit by taking a~=10 
in (3). This value agrees qualitatively with esti- 
mates for a in bulk samples. (In bulk tin, for 
example, the coherence length is 2.5 107° cm 
and the penetration depth 5x10°® cm, so that a 
is approximately eight.) A more quantitative 
comparison with (2) is hardly appropriate since 
its derivation employs the correlation functions 
of bulk material. In a detailed analysis the cor- 
relation functions of the small sample would have 
to be determined and the magnetic field variation 
then derived. (We note that if the boundary of the 
sample could be treated as incoherently scatter- 
ing electrons, the effect of finite size would be 
more important than magnetic field variation in 
diminishing the correction to the free electron 
spin susceptibility. Specifically, if the size, or 
mean free path, 1, were less than 37, R would 
become® 1 - 2//é,. While this treatment of the 
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boundary is known to be unsatisfactory,® and the 
yse of bulk correlation functions is not justified, 
the argument underscores the possible magnitude 
of size-dependent corrections.) 

In the small-attenuation or London limit (a<<1), 
the ratio (2) reduces to the form previously cal- 
culated by Yosida’®: 





exp[ y? + (34)?)”? 
=2 : 4 
R(0, B) J dy {i+expl y?+ (ga)? }?? (4) 
oO 

This limit, which predicts a rapidly decreasing 
Knight shift with decreasing temperature, is ap- 
propriate only very near the critical tempera- 
ture, 1/8,, where it reduces to 


(0, B=) = 1 -0.60[A(g)/A(0) F. (5) 


We are thus led by (2) to predict tentatively the 
following behavior for the spin susceptibility as 
the temperature is decreased: a sharp drop ac- 
cording to (5) immediately below the critical 
temperature, followed, perhaps, by a rise on 
entering the Pippard region, and then a smooth 
decrease to the asymptotic value (3). Observa- 
tions of a similar structure have been reported 
by Knight," but his preliminary data suggest that 
the sharp valley lies immediately above the crit- 
ical temperature. 

We therefore conclude that the Knight shift ex- 
periments are not at variance with the BCS 
theory. Furthermore, although the effects of size 
may be significant, the predictions for bulk sam- 


ples agree qualitatively with the observations on 
small samples. 
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KNIGHT SHIFT IN SUPERCONDUCTORS 


J. R. Schrieffer* 
Institute for the Study of Metals, Chicago, Illinois, 
and John Jay Hopkins Laboratory for Pure and Applied Science, General Atomic, San Diego, California 
(Received September 3, 1959) 


The experimental values of the Knight shift 
found by Reif’ on superconducting mercury col- 
loids and by Androes and Knight? on evaporated 
tin plates when extrapolated to T =0°K show that 
the electronic spin susceptibility x,(0) does not 
vanish for the superconducting ground state. 

These experiments are in disagreement with the 
calculation of Yosida® based on the BCS theory of 
superconductivity* which predicts an exponen- 
tially vanishing susceptibility near T =0°K for a 
uniform magnetic field due to an energy gap 2€,(0) 


for the creation of a pair of quasi-particles with 
total magnetic moment 2y p. Anderson® and 
Rickayzen® have generalized Yosida’s calcula- 
tion to include the effects of a space-varying mag- 
netic field and they find a nonlocal susceptibility 
which is finite at T=0. Heine and Pippard’ have 
suggested an alternative form for the matrix ele- 
ments which enter the theory such that a finite 
Knight shift is obtained for a uniform field; how- 
ever, no one has been able to construct wave 
functions which lead to these matrix elements. 
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A nonvanishing value of x,(0) can be obtained 
by including a mean free path, 7, in constructing 
the normal state wave functions. Let the single- 
particle wave functions in the normal state, in- 


cluding the effects of the scattering centers which 


give rise to the finite mean free path, be y,, 
having energies €,. The Hamiltonian describing 
the interaction of the magnetic field H, which is 
taken to be in the z direction, and the electron 


—— 
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m 
or T 
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where y,, =, and the c’s obey the usual Ferm; 


anticommutation rules. If the BCS states are cop- 


structed in terms of the Yn ’s, the shift in energy 
to second order in H,,,, when averaged over all 





positions of the scattering centers may be ex- 
pressed as 


spins is 





* = rs r . pe pe > ’ 
ai : og es fe, (HEY (Fy )H(F Wp, At )d°r d’r Lle 5 €, ») ay? () 


where L{é, é,” is given by (5.16) of I. In order to perform the sums in (2), we require averages of 
the form 





* i: =, *2. (2 _sinkR sink’R -R/l - < 
‘ic “nt const"™’ Fy, MW, (F WF) ay RRR?” ; aethadtithd, @) 
where the average has been approximated by the 
product of the averages of the m and n’ sums and where 
we have used the results of Weiss and Abrahams® 111 
for the averages involved. The states y, reduce i hg 


to linear combinations of plane-wave states with 
wave number |k/ in the limit /-, 
The nonlocal susceptibility is given by 


For ]~ this expression reduces to the result of 
Yosida; however, for finite 7 one obtains a non- 
vanishing magnetization at T=0. To obtain agree- 


-0°(E™) Av ment with the experiments of Reif and Androes 
X(R) = One) and Knight one would have to choose //£, to be 
0.50 and 0.37, respectively. The coherence 
sinkR sink'R -R/l length, &, for tin is 2.5x10~° cm,‘ so that 
=p 7D e Lie,,€,,). (4) Btn, So» . . 
B kk'R? k’ k 1~107° cm for the latter experiment. This value 


kk’ 
is an order of magnitude larger than the mean 


dimension of Androes and Knight’s plates, which 
suggests that a large amount of specular scatter- 
ing occurs at the surface. The same situation 
exists for Reif’s experiment. 

It is possible that in small superconducting 
samples the low-lying continuum of magnetic 
states which gives rise to the Knight shift can 
be described by a continuous change of pairing 
with magnetization. For example, one could 


Since L(€,, Ep») is nonzero only when |€,|,!€,/! 
Z€,, it is easily seen from (5.3) and (5.15) of I 
that x(R) is related to the kernel K(R) for orbital 
diamagnetism by 


x(R) =x, [5(R) - > ,*(0)K(R)], (5) 


where x, =2u p°N(0) is the spin susceptibility in 
the normal state and K(R) is given by 


JR, Tye -R/l -R/é start with the magnetized state in the normal 
K(R) GAD,’ J(R, T) =e ° (6) metal and begin pairing states on the spin-up 
L " Fermi surface with those on the spin-down Ferm! 


surface in such a manner that a maximum num- 
ber of pair states are connected by the two-body 
interaction, consistent with the conservation 

laws of the problem. If these states being paired 
were strictly plane waves, then momentum con 


Thus, for magnetic fields which vary slowly over 
distances of the order of ~, we have 


g A (0) P 
= P’\d3 yr? - Tr a 
MG) = {xR )dr =) =a | Len 
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servation would forbid the states interacting with 
each other, since in general the pairs would 

have different total momentum. However, the 
experiments mentioned above have been per - 
formed on samples whose dimensions were such 
that the spread in momentum of the single-parti- 
cle states due to finite sample size was large 
compared to the momentum difference Ak due 

to the change in pairing caused by the magnetic 
field. Indeed, Ak ~&,<1 for the magnetic fields 
employed in these experiments, where £, is the 
coherence length. It would appear that pairing 
exactly time-reversed states as one does for 

the ground state, or pairing states which differ 
slightly from these, should not alter the inter- 
action energy strongly and it is possible that one 
would obtain a finite Knight shift in this manner. 
If this description were valid, then the magnetic 
continuum would be analogous to the continuum 

of current carrying states found by pairing 
single-particle plane-wave states (k+qt, -k+qu). 
The vanishing of the nuclear spin relaxation® rate 
as T-0 indicates that the density of these mag- 


netic states would be small just as in the case 
of the current-carrying states and that an energy 
gap would continue to exist for single-particle- 
like excitations. 

The author is indebted to Dr. John Bardeen and 
Dr. G. Rickayzen for several stimulating dis- 
cussions relating to this problem. 
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KNIGHT SHIFT IN SUPERCONDUCTORS 


P. W. Anderson 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received September 9, 1959) 


Ferrell’ has recently explained the observation 
by Reif* and by Androes and Knight? of finite 
Knight shifts in tin and mercury in terms of spin- 
reversing scattering (a spin-orbit effect) at the 
surfaces of the fine particles used. This theory 
might be questioned because it uses concepts — 
coherence length and plane wave states—which 
are well defined only for bulk samples of pure 
superconductors, while the experiments are done 
on fine particles. This Letter demonstrates that 
the spin susceptibility and Knight shift may be 
calculated using the theory of very imperfect 
(“dirty”) superconductors proposed by the author, 
with no further assumptions. The result con- 
firms Ferrell’s physical reasoning. 

Ina very imperfect metal, plane-wave states 
with fixed spin are no longer good one-electron 
functions, so the theory of reference 4 introduces 
hypothetical exact one-electron functions, 


4 


. “Lig ol” IKo)be 5 energy € (1) 


lt is observed that in the absence of magnetic 
centers the Kramers time-reversal degeneracy 


must be present, so that the distinct state 


= + 
Vin “LE o” \ko) ¥ 9’ 


is also an eigenstate of energy €,- Then it is 
shown that the pairing (n, -n) leads to a BCS state 
with essentially the average energy gap of the 
bulk superconductor, explaining the fact that 
such samples as those of Androes and Knight 
have nearly the bulk 7... 

In the presence of spin-orbit scattering, the 
state y,, must not be an eigenstate of the spin; 
in fact it will normally have no average spin 
component in any direction. Thus the spin oper- 
ator S has now only off-diagonal components, 
causing transitions from states n ton’. If En 
and €,,, usually differ by more than the gap, 
clearly the resulting spin susceptibility is not 
much affected by superconductivity; our theory 
merely expresses this physical fact. 

Using time-reversal symmetry and Hermiti- 
city, it may be shown that in terms of scattered 
function fermion operators Cc,» the spin operator 
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electrodynamics of dirty superconductors. 
The susceptibility is to be calculated from 


x = (8? /h*)2) (01S, |m)(m iS, \0)(E,,-Eo)*. (4) 


If we assume a Bogolyubov® ground state defined 
, " " by the “quasi-particle operators” 

=sh - »o ©-€ = Cm Be 
s; 2c, Cee Oniey © ie) a sue *-ve, 
The energy spectrum of the coefficient Sun’ is in 
principle easily determined from a measurement 7. eT 
of the absorptive component of the susceptibility 
in the normal state: 


* 
ba (5) 


(u» and vy are the usual functions of energy « n> 
only the part of S, proportional to a,,ta,,,7 has 


x'"(w)« 27 ole. “<__," w) a \?=f(w7), (3) matrix elements with the ground state; inserting 
cael (5) into (2) we get 
defining the relaxation function f and the spin ; ; 
relaxation time 7. f might be taken proportional S =3h D u vy -uvp)aa Ss ., 
2_2\-1 P s, n n n n a nn 
to (1+w*°r°)™, for instance. n,n’ 
+C.C. + ata terms, (6) 


We actually use (3) in such a way that the aver - 
age over n is taken only over functions with the 


and x f 4) become 
same energy. This is no limitation; there is no we y Seem £6) — 


reason to expect relaxation times to vary over x =48? > IS uv -u yp P(E +E). (1) 
~1°K from the Fermi surface. It is this observa- n,n’ — ee 

tion that one can use an observed relaxation func- Here E,,= (e,,” +€ ar. Clearly this is zero if 
tion to get the energy dependence of the matrix En =€y', but “otherwise not. 


Since u, v, and € are only functions of En» we 
may use (3) to obtain the final answer in terms of 
the energy gap and relaxation function: 


elements that makes the technique work. A 
similar trick should give accurate results for the 





Fcsacoid f ((€ - €’)th™ Jlu(e)v(e’) - u(e’)v(e) P(E +E’) 


lim Jdefde’--- 
€,~0 


(8) 





os |e 





The quadrature could be done numerically; 
rough approximate results are 


X/x, =1-(2€97M),  EgQ<h/t (9) 


face layer, a negative one in the inaccessible 
interior. This will greatly complicate the size 
dependence. 





x./x, = a(€-€')?) /€g",  €y>h/7 
'R. A. Ferrell, Phys. Rev. Letters 3, 262 (1959). 
I am indebted to Professor Ferrell for a preprint of 


this Letter. 
*F. Reif, Phys. Rev. 106, 208 (1957). 


=gh/e7. (10) 


In the case of Knight and Androes’ measure- 





ments on Sn, (9) gives 2€,7/A=2, and using the 
usual energy gap® this means 7 =1.2x107', cor- 
responding to spin reversal once every 10-20 
collisions with the surface: slightly less often 
than Ferrell found. 

It is interesting that a finite Knight shift will 
be observed for almost any size of sample. Once 
the sample is large enough that the above effect 
does not work, the nonlocal nature of the elec- 
tron susceptibility’ will begin to be felt, so that 
a positive Knight shift will be seen in the sur- 
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CYCLOTRON RESONANCE IN METALS* 


J. C. Phillipst 
Department of Physics, University of California, Berkeley, California 
(Received August 19, 1959; revised manuscript received September 8, 1959) 


Well-resolved Azbel’-Kaner resonances (some- 
times showing as many as m =15 subharmonics) 
have been observed in several metals (tin, Al,”»* 
Zn‘) including Cu.° With H normal to the surface 
of the metal, “dielectric anomalies” have been 
observed* even in the extreme anomalous limit. 
The characteristic feature of these effects is that 
they are more striking than predicted by theory.°~® 
The purpose of this note is to elaborate on an 
idea proposed by Chambers’ and Heine® to ex- 
plain the results obtained for H nearly parallel 
to the surface of the metal and to extend this idea 
to give a qualitative explanation of the masses 
observed with H normal to the surface of the 
metal. 

With warped energy surfaces the electronic 
motion in the presence of a magnetic field can be 
discussed, according to Shockley,® in terms of 
“tubes” in k-space. A tube is a region in momen- 
tum space lying between ky and ky +dky, and E 
and E+dé; for degenerate statistics the effective 
tubes all lie at the Fermi surface and can be in- 
dexed by ky. We can associate with each tube a 
drift velocity Yp defined as 


es 1 ('@ 
vplky) =F I Fan, Be 


where the integration is taken over a complete 
period. From the equation of motion, 

dk edr = 

dtc dt’ 
it can be shown that Vp WH. 

Most treatments of this problem have assumed 
that the relaxation time 7 is independent of ky. 
The surface is then regarded as a boundary 
condition. An alternative formulation which helps 
in interpreting line shapes is obtained by making 
the substitution 

wang DO (1) 
+ Ff Bd 
Here y Dy iS the magnitude of the component of 
the drift velocity normal to the surface of the 
metal, 6 is the skin depth, and 8 is a number of 
order unity, similar to the quantity introduced 
ty Pippard’® in his “ineffectiveness” treatment 
of the anomalous skin effect. From (1) it is clear 


that for 6</ only tubes with vp, = 0 will con- 
tribute appreciably to the surface impedance. 

The most surprising feature of the observed 
Azbel’-Kaner resonances is the large number of 
subharmonics that are resolved. These can be 
explained by assuming that the magnetic field 
makes an angle y with the surface, where 
y25/re=2x107 H (for Cu; here r, is the cyclo- 
tron radius, and H ~10°-10* gauss). In this 
case, as Chambers’ and Heine® have indicated, 
Vpn =¥|Vp! and only tubes with Vp ~0 contribute 
to the resonance. In Cu with 7=10~° sec, 6=107° 
cm, v ,=10* cm/sec, and » =10~, the second 
term in (1) will be as large as the first for "'”D 
=107y F: Thus only narrow tubes contribute to 
the resonance, which accounts for the observa- 
tion of unique masses. 

If a particular Fermi surface is centrosym- 
metric, then vp(-ky) =-vp(ky). In particular 
vp(0)=0. The ky =0 tube is also an “extremal 
mass” tube: 


w (ey) =w (0) +0°k,* ere (2) 


and this tube should contribute especially strongly 
to cyclotron resonances. We also have 
Up (ky) =ck,, ven (3) 

In general conditions similar to (2) and (3) will 
not hold simultaneously for tubes other than kz, =0. 
Similar results may obtain, however, for orbits 
about especially symmetric points of Brillouin 
zones such as the centers of the (111) faces in fcc 
lattices. In these cases ky in (2) and (3) should 
be replaced by (ky -ky°), where ky° refers to 
the symmetry point in question. Whether such 
orbits are possible depends on the orientation of 
H relative to the crystal axes, and these special 
orbits may account for the large-angle reso- 
nances observed in tin.’ 

These remarks suggest that the following fea- 
tures should be observed in Cu. For H parallel 
to [100] the resonance associated with the k,, =0 
tube is observed; here (2) and (3) may contribute 
about equally towards a unique mass. As H is 
rotated in the plane of the surface, the ky =0 
tube will continue to resonate; this produces a 
“normal” branch in the plot of m* vs orientation 
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of H in the plane. We expect that 5 will increase 
as H departs from the [100] axis; this would 
ordinarily mean fewer subharmonics would be 
resolved. The fact that more are resolved can 
be explained by assuming that c increases faster 
than 6. If the Fermi surface overlaps the zone 
edge near the center of the (111) face, m*™ of the 
normal branch should have singularities when H 
is near the [110] direction, since then open orbits 
are possible. The normal branch thus “maps 
out” the Fermi surface. In addition to the normal 
branch, several anomalous branches may appear 
associated with the (111) contact areas [k,,° 
=na~*(1, 1, 1)]. 

Quite different line shapes are expected depend- 
ing on whether C =(y7c/5)? or B=6? is larger. 
When more subharmonics are resolved, so that 
C is larger, the masses are more unique, so 
that the asymmetrical lines predicted by Azbel’ 
and Kaner in the ellipsoidal case should be more 
in evidence. In particular the asymmetry should 
be preserved in the higher subharmonics by (3) 
v.hereas (2) leads to a progressively symmetri- 
cal broadening. 

Chambers” has argued that in the extreme_ 
anomalous limit Z should be independent of H 
for H normal to the surface. In this geometry 7 
in (1) should be replaced by 7/(1+iw’rt), where 
(for circular polarization) w’ =w - >. From (1) 
and (3) it can be seen, in agreement with Cham- 
bers’ conclusion, that the decrease in 7 is ex- 
actly compensated by an increase in the number 
of effective carriers. Under these conditions, 
however, Galt et al.* have observed large effects 
in Zn very similar to those obtained under classi- 
cal conditions in Bi.’ It is possible that these 
effects result from higher-order corrections, 
such as the variation of E in the skin depth, which 


are neglected in the ineffectiveness treatment. 
In any case the observed masses should result 
from tubes with YD =0. With wr =10*, 7 replaced 
by 1/w and » =1, we find as in the transverse 
case v £10» 5. Thus the same masses should 
be observed as in the transverse case; this con- 
clusion agrees qualitatively with the results of 
Galt et al.* in Zn. 

I am indebted to Dr. D. N. Langenberg and 
Mr. T. W. Moore for permission to use their 
results in Cu before publication. I have bene- 
fited greatly from conversations with Professor 
C. Kittel, Professor A. F. Kip, Dr. W. Harrison, 
Professor M. H. Cohen, and Professor M. Tink- 
ham. 
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CYCLOTRON RESONANCE IN COPPER* 








D. N. Langenbergt and T. W. Moore 
Department of Physics, University of California, Berkeley, California 


We are reporting preliminary results of a study 
of cyclotron resonance in copper. Observation 
of cyclotron resonance in a sample of natural 
copper consisting of a small number of single 
crystals was first reported by Langenberg, Kip, 
and Rosenblum.’ The present experiments were 
done with an artificially grown single crystal and, 
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because of improved experimental techniques and 
a more satisfactory sample, have yielded signif- 
icantly more detailed results. These are the 
first observations of a resolved cyclotron reso- 
nance in a monovalent metal. 

The experiments were done at 4.2°K and at a 
microwave frequency oi 24 kMc/sec using a con- 
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yentional microwave bridge. One plane surface 
of the coin-shaped sample was cut within two 
degrees of a (110) plane and electropolished. 
This surface formed one end wall of a TE,,, 
cylindrical cavity. The sample and cavity were 
mounted in a magnetic field so that they could be 
rotated together about an axis normal to the sam- 
ple surface, keeping the magnetic field in the 
plane of the sample. Experiments were done 

both with and without magnetic field modulation. 
The results obtained with the two methods were 
consistent. Most of the results presented here 
are based on the field modulation data because 
these data were generally better. 

The experimental results without field modula- 
tion showed a general decrease in sample sur- 
face resistance with increasing magnetic field. 
This decrease was estimated to be approximately 
30% from zero field to the maximum value of 16 
kilo-oersteds. 

Well-resolved resonances were observed for 
all orientations of the magnetic field in the (110) 
plane. Only single masses were observed when 
the magnetic field was parallel to the (100), (110), 
and (111) axes; however, for most arbitrary 
orientations, at least two masses were well re- 
solved. For some orientations, the anisotropy is 
such that the mass continuity is lost when making 
changes in orientation as small as two degrees. 

Cyclotron masses were calculated from the 
period A(1/H) of the derivative maxima, as shown 
in Fig. 1. The masses for H parallel to the prin- 
cipal crystal axes and their behavior with respect 
to orientation are given in Table I; the vaiue of 
m*(111) is in excellent agreement with the value 
estimated by Shoenberg’ from the temperature 
dependence of the de Haas—van Alphen effect. 
Although the observed masses lay between ex- 
treme limits of 0.5m, and 5.0m, (both values 


Table I. Cyclotron masses for copper with H in (110) 
plane parallel to principal crystal axes. 





Behavior of mass 
with respect to 





Axis Mass orientation 

(100) (1.32 + 0. 02)mMe maximum 

(111) (1.30 + 0. 03)me increasing as 
H rotated 
toward (110) 
axis 

(110) (1.12 + 0.02), minimum 
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FIG. 1. The lower panel shows a recorder trace of 
absorption derivative for H in a (110) plane 10° from a 
(100) axis. Jy ¢ is in a (110) plane and 45 degrees 
from H, 55 degrees from the (100) axis. f=24.47 
kMc/sec. No significance should be attached to the 
envelope of the derivative peaks since the lock-in de- 
tector phase varied with magnetic field. The upper 
panel is a plot of reciprocal magnetic fields at the 
absorption derivative maxima against integers; m* 
=e/|wcA(1/H)), where A(1/H) is the slope of the 1/H 
plot. Note the large “phase” discrepancy between the 
two nearly equal masses. 


observed at the same orientation, with H 40 
degrees from a (100) axis), the masses for most 
orientations were between 1.1m, and 1.4m,g. Ex- 
perimental errors in the masses themselves 
were usually less than 2%; the largest uncer- 
tainty in quoting masses for specific orientations 
where the anisotropy was large arises from small 
uncertaintie$ in orientation. 

For finite relaxation times (we estimate w7 ~50 
in the present case), small corrections should 
be applied to these experimental masses because 
absorption derivative maxima rather than absorp- 
tion minima have been used and because the ab- 
sorption minima themselves do not occur at ex- 
actly the cyclotron fields. No attempt has been 
made to apply these corrections since the latter 
one depends upon the “character” of the mass, 
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i.e., whether the associated carriers are ona 
region of the Fermi surface where the mass is 
a relative minimum or maximum, or ellipsoidal 
(Azbel’ and Kaner*), and this is not yet known 
with any certainty. We believe the net mass 
correction to be small in any case, certainly 
less than 5%, since the two corrections tend to 
cancel one another, w7 is large, and there is 
no detectable curvature in the 1/H plots. 

These resonances are perhaps the first in any 
metal for which the relationship with the Fermi 
surface appears fairly clear and from which 
further detailed information about the surface 
may be derived. Pippard,* on the basis of anom- 
alous skin-effect experiments, has proposed a 
model of the Fermi surface of copper which lies 
entirely in the first Brillouin zone and is very 
nearly spherical except for “humps” which con- 
tact the hexagonal zone boundaries around (111) 
directions; Olson and Rodriguez® also conclude 
that the surface contacts the zone boundary, 
from an analysis of magnetoresistance measure- 
ments. The very complexity of our data seems 
to argue against the existence of a single, closed 
convex surface. If Pippard’s model is at least 
qualitatively correct, the following features 
should be apparent in the resonances: for H 
parallel to a (100) axis, closed orbits are avail- 
able to carriers having small momentum along 
the field, and those having no momentum along 
the field should display a resonance, since by 
symmetry they have an external mass. If the 
resonance due to these “equatorial” orbits is 
followed as H is rotated toward a (110) axis, the 
associated masses should change only slowly 
over the nearly spherical part of the Fermi sur- 
face, but should have singularities at the three 
orientations for which the “equatorial” orbits 
just intersect the edges of the contact areas. 

Our data are in excellent agreement with these 
comments. One of the expected mass singulari- 
ties has been observed in detail with H 18°+1° 
from a (110) axis, as compared to 143° predicted 
by the Pippard model; the present data are not 
complete enough for us to observe the other two 
expected singularities. 

For H along the (100) and (110) axes only 8 
or 9 harmonics are resolved and the shape of 
the derivative peaks is fairly symmetrical, while 
for H along arbitrary axes 12 harmonics are 
often observed and the peaks have the asymmetri- 
cal “N” shape characteristic of the Azbel’-Kaner 
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theory for ellipsoidal (unique) masses and high 
wT; in general, the peaks become increasingly 
symmetrical as the order of the harmonic in- 
creases. The implications of these results are 
considered by Phillips.® 

An interesting unexplained feature of the reso- 
nance is that at some orientations the “phase” of 
the oscillations is not that predicted by Azbel’ 
and Kaner; i.e., if the period of the oscillations 
is fitted to the Azbel’-Kaner curve, the positions 
of the peaks are shifted with respect to the theo- 
retical curve. This phase shift is most easily 
observed by noting the intercept of the 1/H plots 
(see Fig. 1). In at least one case a continuous 
and rather rapid variation of the phase is ob- 
served as the magnetic field orientation is changed, 
A phase discrepancy has also been observed by 
Aubrey and Chambers’ in bismuth and by Galt® 
et al. in zinc. 

This work is being continued and a complete 
report will be published at a later date. 

We should like to express our gratitude to 
Dr. J. E. Kunzler of the Bell Telephone Labora- 
tories for providing the copper single crystal 
used in this work and for making a residual re- ) 
sistivity measurement on it, to Dr. J. C. Phillips 
for many helpful discussions and for making his 
manuscript available to us before publication, and 
to Professor A. F. Kip for his guidance during 
the preliminary stages of this work and for his 
critical comments and counsel while the manu- 
script was being prepared. 
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THEORY OF THE RESISTANCE MINIMUM IN DILUTE PARAMAGNETIC ALLOYS 


A. D. Brailsford and A. W. Overhauser 
Scientific Laboratory, Ford Motor Company, Dearborn, Michigan 
(Received September 3, 1959) 


The presence of a minimum in the electrical 
resistivity as a function of temperature, at low 
temperatures, appears to be a general property 
of dilute paramagnetic alloys of a transition 
metal in a noble metal.' Recent experiments’ on 
Cu-Fe alloys indicate that the resistance mini- 
mum is a true bulk property depending on the 
presence of paramagnetic ions in random solid 
solution. 

It has been shown that the scattering of conduc - 
tion electrons by paramagnetic ions yields a 
temperature -independent contribution to the re- 
sistivity.* However, if the Zeeman levels of the 
ions are split by a magnetic field (for example), 
this contribution becomes temperature dependent. 
The anomalous resistivity of paramagnetic alloys 
suggests the presence of scattering centers with 
closely spaced energy states. Such centers occur 
when the interaction between nearest-neighbor 
pairs of paramagnetic ions is considered. If the 
spin of each ion is S, then states of the ion-pair 
are characterized by the total spin J, which has 
integer values between 0 and 2S. The energy 
separation of these states of different J] depends 
on an exchange integral W. We have calculated 
the resistance arising from scattering by such 
ion-pairs, and find that if the interaction of a 
pair is ferromagnetic (W>0), a resistance mini- 
mum should occur. Since this mechanism in- 
volves pairs of ions, the size of the minimum 
should be proportional to the square of the con- 
centration of the transition metal. This quadratic 
dependence appears to be qualitatively confirmed 
for the most dilute Cu-Co alloys investigated by 
Jacobs and Schmitt. [One expects this depend- 
ence to obtain only for extreme dilution, where 
the (2/+ 1)-fold degeneracy of each level is not 
split by other interactions. | 

A mechanism for a temperature-dependent re- 
sistivity is described in the heuristic model of 
Schmitt.° The resistivity caused by elastic scat- 
tering will be temperature dependent since the 
occupancy of the different energy states, having 
different scattering cross sections, varies with 
temperature. Inelastic scattering is temperature 
dependent for the same reason and also because 
the available final states for the scattered elec- 
tron depend on temperature. 


We shall treat the s -d exchange scattering of - 
conduction electrons by ferromagnetically coupled 
pairs. The model Hamiltonian is the sum of two 
terms: 


H,=-WS,°S,, (1) 
and 


H, = V(t) +V(f-R) - 2J(F)8-S, - 2(F-R)S-S,. (2) 


S, and s, are the spin operators for the ions, $ 
the spin operator for an electron in the conduc- 
tion band of the solvent metal, and R is the near- 
est-neighbor distance. V(r) andJ(r) denote the 
spin-independent and spin-dependent coefficients 
of the interaction between a conduction electron 
and an ion. The eigenvalues of H, are given by 


E,=2W((l+ 1) - 288+ 1)], 


each level being (2/ + 1)-fold degenerate. 

We find that the average elastic scattering 
cross section for each of the 2/+1 states of en- 
ergy E, is proportional to J(7+1). If W>0, states 
of lower energy have larger elastic scattering 
cross section, so that this contribution to the 
resistance increases as temperature decreases. 
Therefore a resistance minimum would appear 
to be easily explained. However, one must take 
account of inelastic scattering, which becomes 
frozen out at low temperatures and yields a de- 
creasing contribution to the resistance as tem- 
perature decreases. The net temperature de- 
pendence depends on the close competition of 
elastic and inelastic processes, so one must ex- 
plore the model in detail to show that the elastic 
contribution predominates. 

The scattering arising from the V terms of (2) 
is temperature independent and need not be con- 
sidered further. All cross terms between V and 
J terms may be shown to vanish when square 
matrix elements are averaged over conduction 
electron spin states. The J terms yield the fol- 
lowing elastic and inelastic contributions to the 
resistivity: 


Pol ¥ on Fs (8) 
Pine =*F_8_: (4) 


331 





VoLUME 3, NUMBER 7 


PHYSICAL REVIEW 


LETTERS OcToBER |, 1959 





where 
a= m*N, /aane*h'K®, 


Np being the number of pairs per cc, m the elec- 


tron concentration, and K the Fermi wave number. 


The factors B, are the following integrals of the 
Fourier transform J(q) of J(r): 


2K 
B, -[ J (q) ?[1 + (singR /gR) ]q°dq. (5) 
0 


The second term in the square brackets of (5) 
represents the interference of the scattered 
waves from the two atoms of a pair (averaged 
over random directions of R). One should ob- 
serve that the interference is constructive for 
elastic scattering, whereas it is destructive for 
inelastic scattering. It is precisely this feature 
which allows the elastic scattering to predomi- 
nate over the inelastic. The remaining factors, 
f,, of (3) and (4) contain the temperature depend- 
ence of the resistivities and are 


2S 
f= p (21+ 1) + 1)/4, (6) 
+ J=0 


2s 
f_= 2 byP,_ (2S + 1)?-17\V/(b +P, 4); (7) 


where p, is the equilibrium probability that each 
state of energy E; is occupied. 

Consider now the temperature dependence of 
the resistivity contributions. For W/kT<1, 


i 


Eqs. (6) and (7) become 
f= 4S(S +1)[1+S(S + 1)W/3kT]. 


Therefore, the net temperature-dependent con- 
tribution to the resistivity is 


Ap =aS*(S +1)(B -B_)W/6kT. (8) 


Since B, > B_ [assuming only that J(q) is well 
behaved], the resistivity will increase with de- 
creasing temperature. One can show that this 
increase will be monotonic all the way to 0°K as 
long as 


B /B_>2(4S + 1)/(4S-1). (9) 


If the inequality (9) is not satisfied, the resistiv- 
ity will go through a maximum before reaching 
its 0°K value. In either case the competition be- 
tween (8) and the phonon resistivity should pro- 
duce a minimum. 

The magnitude of the anomalous resistivity in- 
crease also agrees with experiment for reason- 
able values of J and W. A detailed account of 
this work will be published elsewhere. 
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ATTENUATION OF SOUND IN A GERMANIUM CRYSTAL AT ULTRA-HIGH 
FREQUENCIES AND LOW TEMPERATURES* 


E. Roland Dobbs,f Bruce B. Chick, and Rohn Truell 
Metals Research Laboratory, Brown University, Providence, Rhode Island 
(Received June 29, 1959) 


We have measured the ultrasonic attenuation of 
compressional and shear waves in a high-purity 
crystal of germanium at frequencies up to 650 
Mc/sec and at temperatures down to 1.5°K. At 
room temperature there is some evidence for 
the dislocation loss mechanism, but at low tem- 
peratures the attenuation is very small. 

Previous measurements of ultrasonic attenua- 
tion at high frequencies in germanium have been 
at room temperature, and limited to 90 Mc/sec 
in the guided-wave method! and to 300 Mc/sec in 
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the usual, unbounded medium, method.? We have 
added an ultra-high-frequency pulsed oscillator 
to our equipment to extend the frequency range 
to 700 Mc/sec. The low temperatures were a- 
chieved in a conventional liquid helium cryostat, 
the electrical energy passing down a bifilar line 
to the transducer, which was bonded to the par- 
allel-sided specimen. 

The specimen was of n-type germanium, with 
room temperature resistivity 45 ohm-cm and 
net donor concentration 1x10"*/cc. Its flat faces 
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were oriented to be normal, within half a degree, 
to the [100] direction. They were both polished 
toa mirror finish and were parallel within 

95x 10° in. /in. 

The measured (or apparent) attenuation coeffi- 
cient includes, in addition to the absorption in 
the specimen, a loss on reflection at the end 
faces*»* and a loss due to diffraction of the beam.° 
Inthe room temperature measurements we have 
applied a measured reflection correction® and a 
calculated diffraction correction® to both sets of 
results, giving the intrinsic attenuation curves 
of Fig. 1. For the compressional waves, the 
small thermoelastic loss® in the specimen has 
been subtracted, leaving only losses due to other 
intrinsic mechanisms. 

The results for compressional waves are in 
agreement with those of Redwood at the lower 
frequencies and show an approximately square- 
law dependence on frequency—the measured 
slopes being 1.85 for the compressional waves 
and 1.90 for the shear waves. The dislocation- 
damping theory’»? predicts a square-law frequency 
dependence up to a frequency w,,, above which 
the attenuation due to this mechanism should be- 
come nearly independent of frequency. The re- 
sults at the highest frequencies seem to confirm 
the presence of a dislocation resonance, with 
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FIG. 1. Intrinsic attenuation of sound in high-purity 
germanium as a function of frequency at room temper- 
ature . 


Wy, about 500 Mc/sec, corresponding to a lower 
resonant frequency. 

Most of this high-frequency attenuation disap- 
pears at liquid helium temperatures, however, 
as is shown in Fig. 2, for compressional waves 
at 340 and 500 Mc/sec and for shear waves at 
333 Mc/sec. The apparent attenuation is inde- 
pendent of temperature below 15°K and in this 
region is so low that it may be largely due to the 
reflection losses, which we have not yet meas- 
ured at low temperatures. The points shown 
above the dashed curve in Fig. 2 were observed 
on warming the sample, but not on cooling, and 
so may be due to changes in the bond. That a 
satisfactory high-frequency bond was achieved at 
helium temperatures is shown by the photograph 
(Fig. 3) of the compressional waves at 340 
Mc/sec. 

The rise in attenuation that begins at about 20°K 
is similar to that found in quartz by Bommel and 
Dransfeld.* This rise suggests that the phonon- 
phonon scattering associated with Umklapp pro- 
cesses, first discovered in thermal conductivity 
measurements in dielectric crystals, may be 
present here. It is of interest to note that the 
phonon mean free path for Umklapp collisions in 
germanium, as estimated by the method of Ber- 
man et al.,® is about 12 microns at 20°K, or 
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FIG. 2. Apparent attenuation of sound in high-purity 
germanium as a function of temperature. The dashed 
portions indicate regions where experimental effects, 
such as bond transformations, may be causing excess 
attenuation. 
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FIG. 3. Initial pulse and a series of echoes formed 
by 340-Mc/sec compressional waves in a germanium 
crystal at 2°K. The exponential curve is generated 
electronically and used in the measurement of the 
attenuation coefficient. 


roughly the same as the acoustic wavelengths, 
which are 9.8, 14.5, and 10.4 microns for the 
waves shown in Fig. 2. Measurements at still 
higher frequencies are clearly possible at low 





temperatures and we are extending our results 
into the kilomegacycle region. 

We are grateful to Dr. R. Logan of the Bell 
Telephone Laboratories for kindly supplying and 
analyzing the germanium crystal, and to N. 
Pitula and W. E. Oates for orientating and polish- 


ing it. 
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HOT ELECTRONS AND CARRIER MULTIPLICATION IN SILICON AT LOW TEMPERATURE 


W. Kaiser and G. H. Wheatley 
Bell Telephone Laboratories, Murray Hill, New Jersey ) 
(Received June 8, 1959) 


Under equilibrium conditions the mean energy 
of the electrons in a semiconductor is of the 
order of kT,, where T, is the lattice temperature. 
In recent years considerable experimental work 
has shown that the electron energy can be in- 
creased substantially above kT, when an electric 
field is applied to the sample.’ The electrons 
are heated while the temperature of the lattice 
remains constant. For very high electric fields 
the collision frequency between electrons and 
lattice increases and limits the electron drift 
velocity. These effects were first observed by 
Ryder and Shockley in n -type Ge at fields of 
10° v/em.* The increased carrier mobility at 
low temperatures reduces the electric field nec- 
essary for the observation of hot carriers. At 
temperatures where most of the carriers are 
frozen out on impurity levels, new effects are 
observed when the electrons become energetic 
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enough to ionize donor impurities. Carrier mul- 
tiplication by this process has been studied in 
Ge in the temperature range of liquid He.* 
We have measured the electrical resistivity 
and the Hall coefficient as a function of the elec- 
tric field in phosphorus-doped silicon (E ; = 0.044 
ev) at 20°K. The samples (usual size 10x2x2 
mm*) were provided with n* contacts. The car- 
rier concentration, n, and the Hall mobility, 14, 
were calculated from the relations n = 1/Rye and 
kh =Ry/p, where p=E/j. In Figs. 1 and 2 the 
current density j, n, and » are plotted as a func- 
tion of electric field for two samples containing | 








less than 10'* oxygen atoms per cm*. The room 

temperature resistivities and the donor and 

acceptor concentrations are listed on the figures. 
To simplify the discussion we divide our field 

range into four parts: 

(1) At low electric fields, E<5 v/cm, the sam- 
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ple is Ohmic, i.e., the mobility 1, and carrier 
concentration nm, are independent of the applied 
electric field. 

(2) Above ~10 v/cm, the carrier concentration 
increases (approximately by a factor of three in 
Fig. 1). The mobility increases in samples where 
Uy is affected by the scattering of ionized impuri- 
ties (Fig. 2) and decreases slightly where 1, is 
mainly determined by lattice scattering (Fig. 1). 

(3) When the electrons have attained a critical 
velocity* v..=1E =2x10° cm/sec, the mobility 
decreases strongly in all samples investigated. 
The rate of carrier formation tends to decrease. 
This field range is limited by the onset of (4) 

(e.g., in Fig. 2, region 3 is hardly observable). 

(4) At high fields, the carrier concentration 
increases rapidly, resulting in a sharp rise of 
the j-E characteristic. The following observations 
should be added: (a) The breakdown phenomenon 
is a true field effect, as was substantiated on 
samples between 1.1107? cm and 1.2 cm of 
length. (b) The breakdown field E ,, is determined 


by the bulk properties of the specimen. Variation 
(factor 3) of surface-to-volume ratio did not af- 
fect Ep. (c) Ep appeared to be independent of 

the contacts used. Diffused n* layers or alloyed 
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FIG. 2. Current density j, carrier concentration n, 
and Hall mobility » versus electric field. (odc, A 
pulses.) p=0.47 ohm cm. 


(Au-Sb) contacts gave identical results on the 
same material. (d) A specimen was constructed 
which consisted of a large part (length 1 cm, 
cross section 2.9 cm’) with a small needle ex- 
tension on one end (length 0.2 cm, cross section 
2.3x107° cm’). Diffused contacts were provided 
to the bottom end of the large part (A), its junc- 
tion with the needle (B), and to the end of the 
needle (C). Current was passed through contacts 
(A) and (C) while the potential could be measured 
between (A) and (B) or (B) and (C). Breakdown 
was observed in the needle while the large part 
of the crystal was still in the Ohmic range. The 
results were found to be independent of polarity. 
Our observations can be explained as follows. 
Region (1): At low electric fields the electrons 
are in equilibrium with the lattice and the re- 
sistivity is constant. (Note that the contacts are 
Ohmic.) Region (2): The increase in electron 
concentration with electric field is considered 
to be the result of the decreasing thermal recom- 
bination rate (between free electrons and ionized 
donor states) with growing electron velocity. 
Experimental evidence for this effect was re- 
cently given by Koenig for n-type Ge in the tem- 
perature range of liquid helium.*® In samples 
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where ionized impurity scattering dominates the 
carrier mobility, », increases with electric 
field since scattering of ionized impurities is 
less effective for hot electrons.® Region (3): 
Shockley’s theory,’ which considers the inter- 
action of hot electrons with longitudinal acousti- 
cal phonons, predicts a critical drift velocity of 
v,.=1.5xc =1.4x10° cm/sec,*® in agreement with 
the observed value of 2x10° cm/sec. The meas- 
ured mobility decrease does not follow the pre- 
dicted E~”? dependence, however. In n-type 
Ge (T, =300°K) agreement between experiment 
and theory has been reported after consideration 
of interactions with optical phonons.’ No theory 
has been advanced which takes into account the 
multivalley structure of n-type Si. Region (4): 
Several processes have to be considered which 
can give rise to a carrier increase in a semi- 
conductor: Injection of carriers was excluded 
by large area n* contacts. The experiment dis- 
cussed under 4(d) is believed to be conclusive 
evidence that carrier injection through contacts 
is not responsible for the observed breakdown 
phenomena. Field emission from the valence 
band or from the donor levels is ruled out since 
the probability for this process is too small for 
electric fields below 10° v/cm. Heating of the 
sample was avoided by immersing the specimen 
into liquid H, and by using short (square) voltage 
pulses. In the measurements shown in Fig. 1 
and 2, the voltage and current pulses (5 107° 
sec duration; 1 pulse per sec) were found to be 
perfectly flat.'° At higher power level the cur- 
rent increased during the application of the pulse 
due to a small heating of the specimen. For 
pulses of 5x10~° sec duration (15 pulses per sec), 
however, the current remained constant up to 
current densities several decades larger than 
shown in the figures. The process most likely 
responsible for the sharp rise in carrier con- 
centration is the ionization of donor impurities 
by hot electrons.'' Preliminary experiments 
indicate that the breakdown voltage depends upon 
Np and Ny, the concentration of the majority 
and minority impurities. Further investigations 
are necessary to elucidate this point. 

Heating of holes and rapid increase in carrier 
concentration were also observed in boron-doped 
silicon (E; =0.045 ev). The effects are similar 








































to those discussed above; v, was found to be 
1 x10° cm/sec. 

In a recent article, Lampert, Herman, and 
Steele’* pointed out that in silicon the energy of 
the intervalley acoustical phonons is smaller 
than E;, the impurity activation energy of phos- 
phorus. They suggest that low-temperature 
breakdown will not occur, because the hot elec- 
trons lose a large fraction of their energy through 
emission of a single phonon before they reach the 
E; value. The results reported here indicate 
that in spite of this loss mechanism enough elec- 
trons reach the threshold energy £; and carrier 
multiplication can take place.'* For p-type sili- 
con (boron) the optical phonon energy (E opt 
= 0.063 ev) is considerably larger than E ; = 0.045 
ev and impact ionization is expected" in agree- 
ment with our observations. 





'For references see J. B. Gunn, Progress in Semi- 
conductors (John Wiley and Sons, Inc., New York, 
1957), Vol. 2: and S. H. Koenig, J. Phys. Chem. 
Solids 8, 227 (1959). 

*E. J. Ryder and W. Shockley, Phys. Rev. 81, 139 
(1950). 

‘N. Sclar and E. Burstein, J. Phys. Chem. Solids 2, 
1 (1957): S. H. Koenig and G. R. Gunther-Mohr, J. 
Phys. Chem. Solids 2, 268 (1957). 

'The assumption Hall mobility =drift mobility is 
made, which will be correct within 10% in n-type Si 
for magnetic fields of 2 x 10° gauss (C. Herring, pri- 
vate communication). 

°S. H. Koenig, Phys. Rev. 110, 988 (1958). 

SE. M. Conwell, Phys. Rev. 40, 769 (1953). 

'W. Shockley, Bell System Tech. J. 30, 990 (1950). 

8¢=9.2 10° cm/sec is the velocity of the longitudi- 
nal acoustical waves. McSkimin, Bond, Buehler, and 
Teal, Phys. Rev. 83, 1080 (1951). 

*E. Conwell, J. Chem. Phys. Solids 8, 234 (1959). 

4 closer inspection on the microsecond level shows 
some structure at the beginning of the pulse in the 
initial range of the breakdown characteristic. 

'lHot electrons and breakdown resulting from impact 
ionization of donors were recently observed at 1.3°K 
by G. Feher, Phys. Rev. Letters 3, 135 (1959). 

'2Lambert, Herman, and Steele, Phys. Rev. Letters 
2, 394 (1959). 

'Nn the high electric fields of p-7 junctions, electrons 
can gain energies of 1 to 2 ev, far above the energy of 
acoustical and optical phonons. Avalanche breakdown 
by this mechanism has been extensively investigated. ' 
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APPARENT STRUCTURE ON THE FAR INFRARED ENERGY GAP 
IN SUPERCONDUCTING LEAD AND MERCURY* 


D. M. Ginsberg,? P. L. Richards,! and M. Tinkham 
Department of Physics, University of California, Berkeley, California 
(Received September 8, 1959) 


We have performed two different experiments 
which provide evidence for apparent structure 
on the edge of the energy gap in superconducting 
lead and mercury at low temperatures (~1.5°K). 
In one experiment, described in a previous Let- 
ter by two of us,’ we measure the far infrared 
power reaching a bolometer after many reflec- 
tions from the walls of a nonresonant cavity. The 
power is called Pg and Py for the superconduct- 
ing and normal states, respectively. Figure 1 
shows the data for a lead and a mercury cavity. 
The sudden decrease in (Ps -Px)/Py with in- 
creasing photon energy signifies the onset of 
absorption in the superconducting state due to 
the excitation of electrons across the energy gap. 
We believe that the dip in the curve at photon 
energies just below the main absorption edge 
represents a precursor absorption in the super- 
conducting state, since the absorption in the 
normal state is expected to rise smoothly and 
monotonically with increasing frequency. This 
dip is evident in our previously published data’ 
on lead. Repetition of the experiment has con- 

















20 T T qT ' T | q T 
a 

16 4 
a ° LEAD 
3 Ff & MERCURY + 
> 
& 2b 4 
x 
» = 4 
ia) 
© et ‘ 
* 

4b of 
a 

J 

50 


0 10 20 30 40 
FREQUENCY (cm™') 

FIG. 1. Frequency dependence for mercury and lead 
of the fractional difference between the power reaching 
the bolometer in the superconducting and normal states. 
These curves have been normalized for display pur- 
poses so that the ordinate of the lowest frequency point 
is the same for each of them. The horizontal bar 
indicates the approximate bandwidth of the lowest 
frequency point. The bandwidth for the other points 
was ~10%, 


firmed this feature in lead, and has shown an 
even more pronounced dip for mercury. 

The structure was actually first noticed in 
measurements by two of us (D.M.G. and M.T.) 
of far infrared transmission through films about 
10 A thick. These experiments are similar to 
those performed on lead and tin by Glover and 
Tinkham,’ but with accuracy improved by changes 
in optics and radiation detection. The magnitude 
of the ratio Ts/Ty of the transmission through 
the superconducting film to that through the 
normal film is measured as a function of photon 
energy. The data are analyzed, using the Kra- 
mers-Kronig relations and a conductivity sum 
rule,* to obtain the frequency dependence of 
0,(w)/oy and o,(w)/oy, where o, and o, are the 
real and imaginary parts of the superconducting 
conductivity, and o,, is the normal conductivity. 
The simplest models for an energy gap, of width 
hw,, predict that o,(w) is zero for w< Wy, and 
that 0,(w) rises monotonically towards On for 
w>w Actually, there appears to be a prelimi- 
nary hump in 0,(w)/on, for lead before the main 
absorption sets in, as shown in Fig. 2. The curve 
for a second lead film has a qualitatively similar 
frequency dependence. Unfortunately, the mer- 
cury films have not been reproducible enough to 
permit a meaningful calculation of the conductivity 
ratio. However, a dip in the transmission ratio, 
corresponding to a hump in o,(w) /On; is clearly 
evident in the transmission curves of two different 
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FIG. 2. Frequency dependence of o,/ oN for a lead 
film with a resistance 197 ohms per square. The 
bandwidth of each point is ~ 10%. 
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mercury films which we have examined. (The 
transmission through the lead films in the normal 
state disagreed slightly with the value calculated 
from the de resistance. The small ambiguity in 
the interpretation of the data which this introduces 
does not affect the qualitative features of the con- 
ductivity curve.) 

It seems likely that the precursor hump in the 
conductivity curves for the films is due to the 
same cause as the dip in the power ratio curve 
for the bulk samples. For both of the metals 
under discussion, the structure in the film data 
and that in the bulk data occur at approximately 
the same frequency, relative to the frequency at 
which the main onset of absorption sets in. This 
structure may be associated with anisotropy of 
the energy gap. Such anisotropy has been ob- 
served in ultrasonic attenuation measurements‘ 
on tin, and could also account for the observed 
nonexponential electronic specific heat, as has 
been pointed out, for instance, by Cooper.® The 
electronic specific heat which is inferred for 
lead and mercury from critical field data shows 
a deviation from the exponential which is an order 
of magnitude greater than that for any of the 
other superconductors measured.*®»? This large 
deviation may well be associated with the marked 
structure that we observe for these two metals. 
On the other hand, the structure may be due to 
the production of collective excitations by the 
absorption of photons. Anderson has pointed out,® 
on the basis of microscopic considerations, that 
collective excited states might have energies 
which lie below the top of the gap predicted by 


the theory of Bardeen, Cooper, and Schrieffer,® 
The collective excitations may perhaps also be 
described macroscopically in terms of longitudi- 
nal oscillations in the metal. In the case of bulk 
samples, transverse oscillations may also play 
a role in explaining the observed effects. Both 
of these possibilities have been proposed tenta- 
tively by Ferrell.’° 

The far-infrared experiments on these and 
other metals will be reported more fully in forth- 
coming publications. 
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OVERHAUSER EFFECT IN METALLIC LITHIUM 


M. Gueron and Ch. Ryter 
Centre D’Etudes Nucleaires de Saclay, Gif sur Yvette, Seine et Oise, France 
(Received August 24, 1959) 


The enhanced nuclear polarization produced by 
the saturation of electron spin resonance in met- 
als can be detected in two ways: (a) by observing 
the enhancement of the nuclear signal, a method 
widely used since the pioneer work of Carver and 
Slichter’; (b) by observing the shift of the elec- 
tron spin resonance brought about by the polari- 
zation of the nuclei, a method proposed by Over- 
hauser’ and analyzed by Kaplan® but never car- 
ried out experimentally so far. 

We present a preliminary report of a detection 
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of Overhauser nuclear polarization in lithium 
metal using both methods. The sample was lith- 
ium hydride heavily irradiated by pile neutrons. 
A sharp (0.3 gauss) electron spin resonance line, 
observed by Doyle et al.* in LiH irradiated with 
ultraviolet radiation, was attributed by them to 
conduction electrons in colloidal particles of lith- 
ium metal. In our sample, with negligible in- 
homogeneous broadening, the electron spin re- 
sonance line had at room temperature and 10 000 
Mc/sec a width of 0.125 gauss (peak to peak of 
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the absorption derivative). The metallic char- 
acter of lithium was confirmed by the observa- 
tion of the nuclear signal of Li’: the line, with 
a width of 0.4 gauss, is clearly diffusion-nar- 
rowed and presents a relative frequency shift of 
0.024%, equal to the known Knight shift in lith- 
jum metal. The narrowness and the symmetrical 
shape of the electron spin resonance line at 
300°K, exhibiting the unusual purity and small- 
ness of metallic particles in irradiated LiH, 
make this sample particularly suitable for Over- 
hauser experiments in high fields and at low 
temperatures. 

(A) Experiments at 300°K: A sample of LiH 
surrounded by a nuclear magnetic resonance 
coil was placed inside a microwave cavity fed by 
a magnetron delivering up to two watts of conti- 
nuous microwave power. The nuclear magnetic 
resonance of Li’ was observed in a field of 3300 
gauss while the electron resonance was being 
partially saturated. Because of the shielding of 
the sample from the microwave field by the nu- 
clear magnetic resonance coil, no reliable esti- 
mate of the field at the sample could be made. 
The enhancement A of the nuclear polarization is 
given by 

Y i Z, 
: 7-5 =1690=*s, (1) 

YLi? “le le 


A-l= 





where 7, is the nuclear spin lattice relaxation 
time, T,, the relaxation time caused by electron- 
nucleus coupling only and s the saturation param- 
eter of the electron resonance, related to the 
microwave power P by 1/s =1+a/P, where a is 
a constant depending on the apparatus and the 
sample. The maximum signal enhancement ob- 
served was approximately 100, the uncertainty 
being caused by the weakness of the unenhanced 
signal. A plot of 1/A obs 28 4 function of 1/P 
was, within experimental error, a straight line 
from which a value of A > 150 could be extra- 
polated. This would require T,/T,. 70.1 rather 
smaller than the value 7,/T,¢~*0.7 at room tem- 
perature and 3300 gauss deduced from the re- 
sults of Holcomb and Norberg.® This discrepancy 
may be due to the presence of grains appreciably 
larger than the electronic (but not the nuclear) 
skin depth, which would give a full contribution 
to the natural signal but not to the enhanced sig- 
nal, 

(B) Experiments at 4.2°K: The relative shift D 
of the electron spin resonance frequency is re- 
lated to the enhanced nuclear polarization (I,) 


and to the Knight shift K by 
N 1 
Doky (1) i,’ (2) 


where x,» is the electron susceptibility, N the 
number of nuclear spins per unit volume, and 
(1,) is given by 
nas Ay, y fll + 1)H, 
- 3kT ' 


where A is the enhancement (1). With H,=3300 
gauss, T=4.2°K, and Xp =2-1 x107®,® the absolute 
shift expressed in gauss is 


4 =DH,=7.2 Ay,/Ve" (4) 


(3) 


The electron spin resonance was observed at 
10000 Mc/sec for various microwave levels up 
to 7 milliwatts, corresponding to an estimated 
rotating microwave amplitude H,~>75 milligauss, 
and to a saturation parameter s>0.3 for grains 
smaller than the skin depth 6. For grains of 
dimension d>6, and for such a narrow line, it 
can be expected, because of electron spin diffu- 
sion, that the saturation parameter will be uni- 
form within the grain but reduced by a factor of 
the order of (5/d)*,’ resulting in a spread in the 
observed electron spin resonance shifts. As the 
microwave power is increased, two striking ef- 
fects occur (Fig. 1): (a) The width and the shape 
of the line change considerably. (b) The shape 
depends on the direction of the sweep. 

The spread of the curve 1(b), of the order of 2 
gauss, agrees with the value of the shift pre- 





26 ah P=.7 mW anh eas 2. 
V 
3a P= .O7mW WN 3» 
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INCREASING FIELD DECREASING FIELD 


FIG. 1. Recorded electron spin resonance absorp- 
tion derivative in neutron-irradiated LiH for various 
power levels P at 4.2°K and 10000 Mc/sec. One divi- 
sion on the horizontal axis corresponds to one gauss; 
the sweep is 0.4 gauss/min. 
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dicted by Eq. (4) for a saturation parameter 
s *0.3, that is for the smallest grains. 

The dependence of the line shape on the direc- 
tion of the sweep is in qualitative agreement 
with the prediction of Kaplan* if the finite nuclear 
T, (10 sec at 4.2°K) and a spread in the size of 
the grains are taken into account. 

The existence of a distribution of electron spin 
resonance shifts was also exhibited directly: 
after saturating the electron spin resonance line 
for several T,, a large field modulation (5 gauss) 
was introduced suddenly. The electron line being 
then saturated during a small fraction of the 
modulation cycle only, the enhanced nuclear 
polarization could relax back to approximately 
its normal value. The electron line observed on 
the scope presented immediately after the intro- 





duction of the modulation, an asymmetrical 
broadening of the order of 0.4 gauss, which col- 
lapsed into a narrow line with a time constant of 
the order of T,>10 sec. 

The interest of Dr. A. Abragam in this work is 
gratefully acknowledged. 
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ELECTRON COLLISION FREQUENCIES IN NITROGEN AND IN THE LOWER IONOSPHERE* 


A. V. Phelps and J. L. Pack 
Westinghouse Research Laboratories, Pittsburgh, Pennsylvania 
(Received July 6, 1959) 


Attempts’»” to correlate laboratory measure- 
ments of electron collision frequencies in nitro- 
gen with the results of measurements using 
rockets passing through the D-layer of the iono- 
sphere have not been very successful. Our meas- 
urements of the electron collision frequencies 
for thermal electrons in nitrogen, combined with 
an improved analysis of the data obtained with 
the rockets, appear to remove the discrepancies. 

An improved version of the electron drift vel- 
ocity tube used by Bradbury and Nielsen® has 
been used to measure the mobility of electrons 
in nitrogen at such low electric fields that the 
electrons are in thermal equilibrium with the 
gas.* The measured values of the product of 
electron mobility, , and gas density, N, are 
uN =1.10x10%, 3.5107, and 2.8x10" cm™ 
volt~* second~! at 77°K, 300°K, and 373°K, re- 
spectively. Following Phelps, Fundingsland, and 
Brown’ (PFB) these results are analyzed by ex- 
pressing the reciprocal of the momentum trans- 
fer collision frequency, Vm), as a power series 
in the electron energy, u, and substituting into 
the standard expressions for electron mobility® 
to obtain a power series in the most probable 
electron energy, kT/e. The power series for 
the mobility is then fitted to the experimental 
data to obtain the coefficients of the series for 
Vy(u). The resulting v,,(u) is shown by the lower 
solid curve in Fig. 1. The dashed curves show 
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FIG. 1. Momentum transfer collision frequencies 
for monoenergetic electrons in nitrogen as a function 
of electron energy. The energy range of the thermal 
equilibrium experiments is shown by the values of 
kT/e at which the measurements were made and 
serves to indicate the range of validity of the curves 
of collision frequency as a function of energy. The 
mobility data are found to fit the relation pN = 9.6 
x 10% (kT/e)~! - 1.39 x 10" (kT /e)~* em=! volt~ sec” 
within the experimental error so that N/ vyy(«) =8.2 
x 10° (u)-!- 5.9 x 109 (u)-? cm*/sec. 
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the results of analysis of (a) thermal equilibrium 
measurements of microwave conductivity by PFB 
and (b) drift velocity and diffusion measurements 
at average electron energies above thermal by 
Crompton and Sutton.’ The upper solid curve 
shows the results of Anderson and Goldstein.*® 
Figure 1 shows very good agreement between 

the present results and those of PFB and satis- 
factory agreement with the data from Crompton 
and Sutton. These data show that the electron 
collision frequency in nitrogen is very nearly 
directly proportional to the electron energy. We 
shall use this approximation in the following cal- 
culations. 

In order to evaluate the electron collision fre- 
quency in air, we must estimate the collision 
frequency for electrons in oxygen. The results 
of analysis of drift velocity and diffusion meas- 
urements in oxygen show that the collision fre- 
quency for electrons in oxygen is about two-thirds 
of that for nitrogen for energies down to 0.2 elec- 
tron volt.’ We shall assume that the ratio of 
collision frequencies for oxygen and nitrogen 
continues to be significantly less than unity down 
to thermal energies (0.026 ev) so that the error 
is less than 20% when Vn 4) for air is taken 
equal to that for the nitrogen alone. 

The application of the above results to the anal- 
ysis of the properties of the ionosphere is illus- 
trated by the following re-evaluation of recent 
data obtained from rockets fired during a polar 
blackout. First, one notes that because of the 
energy dependence of the electron collision fre- 
quency in nitrogen, relations which assume that 
the collision frequency is independent of energy 
should not be used to analyze measurements 
made in air. In order to facilitate calculations, 
the function & used by Kane’ as a measure of 
the ratio of the differential absorption to the 
refractive index is rewritten in terms of func- 
tions which properly average Vv», (4) over the 
electron velocity distribution.*'° Thus, 


Q(y, Y) = [o (r5) - 00 JIle,0) + 1.40.(y9)}", 


where y=v (u=kT/e)/w, r¢=7/(1-Y), 0 
=7/(1+Y), 0, and o; are the real and imaginary 
parts of the conductivity of the electrons in the 
gas, v(u=kT/e) is the collision frequency for 
monoenergetic electrons at the energy u =kT/e, 
and Y is the ratio of the angular gyrofrequency 
due to the earth’s magnetic field to the angular 
frequency of the radio wave. Plots of Q vs } 

for the value of Y appropriate to Kane’s data are 
shown in Fig. 2 for the cases of collision fre- 
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FIG. 2. The function Q(y, Y) vs y for electron colli- 
sion frequencies directly proportional to energy and 
independent of energy. The curves are calculated for 
Y =0.208 which is the value appropriate to the earth’ s 
magnetic field at the site of the rocket measurements 
and to the frequency (7.75 megacycles/ second) used. 
Values of Q for y>1 are not shown since, for y?>>1, 
Q(y, Y) becomes extremely sensitive to the exact form 
of vy,(u) at small u. 
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FIG. 3. Electron collision frequencies for mono- 
energetic electrons as a function of height above sea 
level. The error limits indicated are calculated from 
those given by Kane and would be approximately the 
same magnitude for the lower set of points. Note that 
since the mean values of Q computed from Kane’ s data 
for heights of 61.3 and 63.1 kilometers lie above the 
curve of Q(y) for Vy,(4) =U, only the lower limits lead 
to real values of v,,(u =kT/e). 
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quency directly proportional to energy and inde- 
pendent of energy. The values of v, (u =kT/e) 
inferred from a given observed value of Q for 
the two different energy dependences of v_ (u) 
are in the ratio of 2.5:1 at very small values of 
v../w. Note, however, that since this ratio 
changes with v,/# one cannot use the Q curve 
for constant v,,, shifted to the left by a factor of 
2.5. 

The points of Fig. 3 show the collision fre- 
quencies for monoenergetic electrons in the D- 
layer computed from Fig. 2 and from Kane’s 
values of Q. The curves give the collision fre- 
quencies for monoenergetic electrons in nitrogen 
at the gas densities and temperatures given by 
Nicolet® using our Vy) (lower curve) and using 
Nicolet’s extrapolation of results of Anderson 
and Goldstein* (upper curve). The agreement 
between the lower curve and the points is within 
the scatter of the rocket data. This agreement 
can be considered as evidence for the correct- 
ness of the gas density and temperature data ob- 
tained from other rocket studies.” 





"This work was supported in part by the Advanced 
Research Projects Agency, the Office of Naval Re- 
search, and the Air Force Special Weapons Center. 
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MOMENTS OF INERTIA OF EVEN-EVEN RARE EARTHS* 


J. J. Griffin and M. Rich 
Los Alamos Scientific Laboratory, Los Alamos, New Mexico 
(Received September 8, 1959) 


Belyaev' has considered the effects of nuclear 
pairing interactions in analogy with the theory of 
superconductivity.?~* His treatment leads to the 
following formula’ for the moment of inertia of 
an even-even nucleus: 


IC kIs_IL) I? 
x f 
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a’+(é, - A)(E_- A) 
Jo 
as a 


5 =h?)> 
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EB, =e, -ar +a’, (2) 


and the sums, / and k, run over all single parti- 
cle states of a self-consistent spheroidal well. 
The “chemical poteritial” is the solution of the 
equation 


DA -(E, -aP a7 [z,-alen. 8) 


N here is the number of neutrons or protons out- 
side of closed shells. The quantity A is equal to 
one-half the induced energy gap. 

Numerical computations based on Eq. (1) have 
been made for twenty-six rare earth nuclei which 
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Table I. Rotational energies of even-even nuclei. Columns one to three list the nuclei considered, their defor- 
mations,2 and their assigned energy gap parameters” which yield the calculated energy of the first rotational (2+) 
state in kev (Column 5). Columns 4 and 6 give the observed rotational energy: © and the percent difference be- 
tween experiment and theory. The energy gap is given in units of w,", which is approximately 41A~'? Mev for a 
nucleus of mass number A. 
















































Eexp FE theor % 
Nucleus 5coulomb A/hw,° (kev) (kev) error 
goNdgg'*” 0.24 0. 139 130 139.5 + 6.8 
gg5IMgq!*” 0.27 0.139 122 123.0 + 0.8 
Sm,,!*4 0.31 0.128 83 99.7 +16.7 
gt Gedy!" 0.28 0.139 123 120.0 - 2.5 
Gdy,'* 0.39 0.128 te) 85.6 - 4.0 
Gd,,'** 0.44 0.117 79 68.3 -15.7 
Gdag,'*° 0.45 0. 108 76 66.1 -15.0 
eeD You 0.33 0.117 86 92.1 + 6.6 
Dygg!*? 0.34 0.108 82 84.4 + 2.8 
Dyyg'*! 0.39 0.100 73 70.0 - 4.3 
esE Teg 0.31 0.108 90 88.3 - 1.9 
Erg,'* 0.100 80 85.9 + 6.9 
Erjoo!®*® 0.314 0.093 80 83.7 + 4.4 
Eres’ 0. 088 79 77.7 - 1.7 
7e¥Di99°** 0.28 0.093 84 86.1 + 2.4 
Ybj99!"? 0. 088 78 80.1 + 2.6 
Yhiee’™ 0.294 0. 085 76 77.1 + 1.4 
Ydje6'™ 0. 085 82 74.7 - 9.8 
r2Hf9,'™ 0.28 0.085 89 83.8 - 6.2 
Hf,.,'* 0.29 0. 085 91 78.5 -15.9 
Hf,e'”* 0.26 0. 086 93 90.3 + 3.0 
uW 108 0.25 0. 086 100 101.9 + 1.9 
Wiso!®* 0.23 0.090 112 110.9 - 1.0 
W312! 0.23 0.096 124 118.3 - 4.8 
16084 39'" 0.19 0.090 137 130.9 - 4.8 
760841! 0.17 0. 096 155 166.8 + 7.1 





See reference 8. 
e reference 7. 
“See reference 9. 
These deformations were obtained from measurements on unseparated isotopes. 


exhibit rotational spectra. The eigenvalues, €,,, Atomic Energy Commission. 
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2L. N. Cooper, Phys. Rev. 104, 1189 (1956). 
3Bardeen, Cooper, and Schrieffer, Phys. Rev. 108, 
1175 (1957). 


certain modifications indicated by a recent study 
of the spectra of odd-mass nuclides.® The energy 


gaps were obtained from a smooth curve drawn ‘N. M. Bogoliubov, J. Exptl. Theoret. Phys. U.S.S.R. 
through the experimental pairing energies of 34, 58 and 73 (1958) [translation: Soviet Phys. JETP 
Johnson and Bhanot,’ and the quadrupole defor - 34(7), 41 and 51 (1958)). 
mations, from Coulomb excitation studies.*® S. G. Nilsson, Kgl. Danske Videnskab. Selskab, 

The results are summarized in Table I. They Mat. -fys. Medd. 31, No. 11 (1959). 


6§. G. Nilsson and B. R. Mottelson, Kgl. Danske 
Videnskab. Selskab, Mat.-fys. Medd. 1, No. 8 (1959). 
"W. H. Johnson and V. B. Bhanot, Phys. Rev. 107, 


indicate that Eq. (1), together with the Nilsson 
model, provides a remarkably good description 


of the moments of inertia in the rare earths, 1669 (1957). 
thereby lending support to the superconductor *K. Alder et al., Revs. Modern Phys. 28, 432 (1956) 
theory of nuclear pairing effects. and references cited therein. 








: *Mihelich, Harmatz, and Handley, Phys. Rev. 108, 
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NEW MODEL FOR VIBRATIONAL SPECTRA IN EVEN-EVEN NUCLEI*? 


T. Tamura! 
University of California, Los Angeles, California 


and 


L. G. Komai 
Hughes Microwave Laboratory, Culver Sity, California 
(Received September 14, 1959) 


The vibrating surface model,’ with its predic- wave functions of each state in the WJ model, it 
tion of equal level spacings and the spin sequence* is not hard to see that the desired change in the 
0, 2, andatriad 0, 2, and 4, is the simplest spectrum will be obtained by introducing a poten- 
among those presented in explaining the vibra- tial energy term of the form® 
tional spectra of even-even nuclei. Experimen- 
tally observed spectra, however, indicate devia- u nin? “B 3 cos"3), (1) 
tions from it among which the more important 
ones are*®: (i) The degeneracy of the triad is + 
lifted and usually’ £.(2)<£,(4); (i) the ratio —" a38___.--2 —______ 3 


E.,(2)/E,(2) is generally greater than two. 

The first approach to explain these deviations 
was proposed by Scharff-Goldhaber and Weneser*® 
who considered the coupling of single particles 
to the simple surface vibrational motions. This 
model gives rule (ii) above, but for (i) it pre- 
dicts E,(2)>E,(4). Recently Raz® found, by ex- = ae = Q—_____—_ 240 
tending the GW calculation, that it was possible 128 ¢ 230 ~-—-~- £ 2.24 
to satisfy also (i); however, as will be discussed i es: nile z = 
below, this treatment appears to possess other 
difficulties. 

Wilets and Jean® proposed a model very similar 
to the simple surface vibrational model, but 
which differs in that a finite, shape-independent 
equilibrium deformation is assumed; (the y- 
unstable model). It explains (ii) but gives E,(2) 
=E, (4). sso ot 100 e. 

Figure 1(a) shows the Cd'"* spectrum obtained 
experimentally by Motz’ in which all members 
of the triad as well as the 3, state were observed; 
(this latter may be considered as a member of 
the third excited surface vibrational state). This 
spectrum bears out (7) and (i) and in addition 
indicates that (iii) E,(0)=E,(4)=E,(2), and (iv) 
the ratio of £,(3)/E,(2) is close to, but some- kev I* 
what larger than three. The calculation of Raz (a) (b) 
predicts too large a value for E,(0) and for E£,(3) 
and thus fails to reproduce the new rules (ii2) FIG. 1. Spectrum of the lower states of Cd'"'. 
and (iv). (a) Experimental; (b) theoretical. On the left-hand 

The WJ model is consistent with (ii) and (iv), side of (a) are given the excitation energy in kev and 


. . . , ae ‘ 5 the spin of each state, while on the right-hand side 
but is inconsistent with (i) and (iii). The incon the ratio E,(/)/E,(2) is given. ‘The 2+ state at 1368 











o g——n00 g— ———_ an 


sistency, however, is not too drastic and this kev seems to be due to a single-particle excitation and 
fact suggests that with a slight modification of the absence of the corresponding state in (b) is without 
the model these discrepancies might well be significance since we here consider only collective 
removed. By examining the explicit form of the states. 
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to the WJ potential 
5C(8 = By)’. (2) 


The spectrum obtained with the appropriate choice 
of the parameters in (1) and (2) is shown in Fig. 
i(b) and it is seen that, in fact, good agreement 
with experiment is found.* The relative intensi- 
ties of electromagnetic transitions between sev- 
eral states are also calculated and compared with 
experiment’? in Table I, and again good agree- 
ment is obtained. 

We now discuss the justification for introducing 
the term (1). For this purpose we first note that 
the WJ potential (2) can be described by curve (a) 
of Fig. 2, where the right and left halves, re- 
spectively, correspond to 7 =0° and > =60° (i.e., 
prolate and oblate deformations). Since the 
potential (2) is )-independent, curve (a) has a 
right-left symmetry. A potential energy curve 
which is consistent with the parametric values 
(summarized in Table II) used in the present 
calculation is shown as curve (b) in Fig. 2. This 
curve has a right-left asymmetry which is clear- 
ly the result of the addition of the term (1) to the 
WJ potential. If we were to derive this type of 
potential from “first” principles, we should have 
taken, e.g., the Nilsson scheme” and added the 
single-particle energies for various values of 
the nuclear deformations. With such a procedure 


Table I. Comparison of the experimental and the 
theoretical values of the relative reduced transition 
probabilities of y rays between several states. All the 
y rays involved are E2 and therefore are not specified. 











Experimental Theoretical 
T(2.— 2;)/ T(2,— 04) 1.37+0.31 1.70 
T(2.— 0;)/ T(2,— 04) 0.020 + 0.003 0.014 





Table II. The parametric values used to obtain the 
theoretical spectrum of Fig. 1(b). In determining 
these values use has been made of the half-life of the 
2; state (see reference 10) as well as the energy E,(2) 
shown in Fig. 1(a). For the m=2 term of (1), itis 
found that a very small value for k, yields good agree- 
ment with experiment; therefore, we show here the 
extreme case of k,=0. 
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FIG. 2. Potential energy curve in various approxi- 
mations. (a) Schematic representation of the WJ 
potential (2). (b) Potential energy curve which is con- 
sistent with the parametric values given in Table II, 
and which produces the theoretical spectrum (b) of 
Fig. 1. 


an asymmetric energy curve is expected and thus 
would seem more natural than a symmetric one. 
To get a curve which can be compared directly 
with curve (b) of Fig. 2, however, it is necessary 
to take into account the residual interaction and 
such a calculation is now under way. 

Regarding vibrational nuclei other than Cd", 
we first note that the ratio E,(2)/E,(2) is less 
than two in several Hg and Pt isotopes. The re- 
cent theory of Davydov and Filippov’” predicts 
values for this ratio that always exceed two, but 
in our model ratios less than two are easily ob- 
tained by allowing negative values for £, in (2). 
The DF theory further predicts that £,(2) +£,(2) 
=E,(3), while in our theory £,(2)+£,(2) is greater 
than or less than E,(3) for E,(2)/E,(2) greater 
than or less than 2, respectively. To the extent 
that only for Cd*** and Pt'* are data available 
[in the region where E,(2)/E,(2) does not deviate 
very much from two], we find that this latter 
prediction agrees with experiment. 

Finally, it should be mentioned that the DF 
theory can say nothing about the 0, state, and 
while £,(2) is correctly given, it predicts a 
value for E,(4) in Cd*"* which is somewhat too 
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large. DF stress the agreement of their theory 
with experiment regarding the relative intensi- 
ties of the M1 and £2 transitions between the 
2, and 2., states, in addition to the relative in- 
tensities of other £2 transitions. However, it 
might be much more reasonable to consider the 
effect of single-particle excitations which begin 
to assume importance at this point.’* This mat- 
ter, however, will be discussed elsewhere to- 
gether with a more detailed account of the con- 
text of the present note. 

The authors wish to thank Professor S. A. 
Moszkowski for his many stimulating discussions. 
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ANGULAR DISTRIBUTIONS FROM STRIPPING REACTIONS OF LOW Q VALUES* 


J. P. F. Sellschop 
Nuclear Physics Research Unit, University of the Witwatersrand, Johannesburgh, South Africa 
(Received June 22, 1959) 


Proton groups from the reactions Li’(d, p)Li® 
and C'*(d,p)C'** (3.09-Mev state) for a range of 
deuteron energies between 0.5 and 2.5 Mev have 
been studied using high-resolution magnetic anal- 
ysis. In fitting the measured angular distribu- 
tions for the protons to a simple stripping theory, 
uncorrected for Coulomb and other perturbing 
effects, an unusually high degree of agreement 
was found. 

Wilkinson’ suggested that such agreement might 
be expected for reactions of fairly low Q values. 

It is well known that the deviations of measured 
stripping (deuteron- induced) reactions from the 
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simple Butler? and Born-approximation*® theories 
for moderate energies (> 10 Mev, say) take a 
characteristic form. In general the first maxi- 
mum can be unambiguously matched, but with 
increasing angle the agreement soon deteriorates. 
In particular the well-defined first minimum in 
the theoretical curves is filled in and subsequent 
maxima are poorly defined. In fact some 20 or 
so degrees after the first maximum the meas- 
ured intensity is normally persistently higher 
than the theoretical prediction. This is attributed 
to competing compound-nucleus formation. The 
perturbing effects of nuclear and Coulomb inter- 
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actions also play their role in distorting the sim- 
ple comparison. om 
At low deuteron energies the deviation of meas- 
ured distributions from simple theoretical ones 7 (dem | 
becomes more serious. This is not unexpected oueeat 


as the effect of compound-nucleus formation and 
interaction effects are anticipated to be more 
prominent. : 

In order to improve agreement between experi- 
ment and theory, many modifications have been 
made to the simpler formalisms. Among these 
can be mentioned the work of Friedman and 
Tobocman,* Daitch and French,°* Horowitz and 
Messiah,® Grant,” Tobocman,® and Tobocman and 
Kalos.° These refinements show that the effect 
of considering the Coulomb interaction is to dis- 
place the angular distributions toward larger 
angles, to broaden the peaks and fill in the val- 
leys, and to reduce the total cross section. On 
the other hand, the effect of introducing nuclear 
interactions is to displace peaks towards smaller 
angles, to sharpen these peaks, and to reduce 
the total cross section. The distortion of simple 
stripping curves which can be achieved is strik- 
ingly great. It is unfortunate that when allowing 
for all possible effects, agreement with experi- 
ment is readily obtained, so that it is impossible 
to say whether the parameters employed to a- 
chieve such agreement have any particular physi- 
cal significance. The conclusion is drawn that in 
order to avoid having to make hypotheses on the 
values of these parameters, experiments should 
be performed at energies well above the Coulomb 
barrier. 

Two examples of the results obtained for the 
Li'(d,p)Li® (ground state) reaction with Q value 
of -0.188 Mev are given in Fig. 1. These were 
obtained for incident deuteron energies of 1.9 
and 1.5 Mev, respectively. 

Another reaction of small Q value has been 
studied, viz., C’*(d,p)C*** (3.09-Mev state). The 
Q value is -0.367 Mev. An example of the re- 
sults obtained in an angular distribution, for an 
incident deuteron energy of 2.1 Mev, is given in 
Fig. 2. 

The theoretical curves associated with the 
measurements given in these diagrams were ob- 
tained from the non-Coulomb, noninteraction 
stripping theory formalism of Friedman and 
Tobocman,* using for convenience the tables and 
graphs of Enge and Graue,*° based on this for- 
malism. Other stripping derivations have also 
been tried, also with striking agreement. 

Excitation functions have been measured for 
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FIG. 1. Angular distribution of protons from the 
reaction Li'(d,p)Li®. (a) Eg=1.9 Mev. The points re- 
present experimental observations; the solid curve re- 
presents a calculated stripping distribution with R=6.0 
x 107") em and /=1. (b) Eg=1.5 Mev. The points 
represent experimental observations; the solid curve 
represents a calculated stripping distribution with 
R=6.5%10-"8 cm, and /=1. 
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FIG. 2. Angular distribution of protons from the 
reaction C!%(d,p)c'* (3. 09-Mev state) Eq=2.1 Mev. 
The points represent experimental observations: the 
solid curve represents a calculated stripping distribu- 
tion with R=7.410-'? em. and /=0. 


both reactions, Li’(d,p)Li® (ground state) and 
C'?(d, p)C'5* (3.09 Mev). In the lithium case one 
resonance was closely studied. Angular distri- 
butions were measured at the resonance energy, 
and at energies closely flanking the resonance. 
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No change-—other than the increase in intensity 
associated with the resonance—in the marked 
degree of agreement between experimental and 
theoretical distributions was found. 

In the carbon experiment a number of resonances 
were found and the effect of a particularly strong 
one in the angular distributions was studied in the 
same way as in the lithium case. Here too there 
was no significant change in the degree of good- 
ness of fit with the simple stripping formalism 
used. 

On the basis of these two experiments it would 
indeed appear that for small-@ reactions at low 
bombarding energies there is particularly little 
influence on the shape of the angular distributions 
from Coulomb and nuclear effects. 

The physical picture underlying this proposal 
is as follows." The momentum of the outgoing 
proton in a (d,p) reaction is made up of about 
half the momentum of the ingoing deuteron and of 
a part of the internal momentum of the deuteron 
at the instant at which stripping occurs. In the 
case of a large-@ reaction the total proton mo- 
mentum must be high. If, however, the proton’s 
share of the momentum of the ingoing deuteron is 
small—when the deuteron energy is small-—the 
large deficit in momentum must be obtained from 
the ground-state wave function of the deuteron. 

In order to get this large latter contribution, the 
separation of proton and neutron in the deuteron 
must be small at the instant of stripping. Fora 
small-Q@ reaction at low energy, on the other 
hand, the contribution to the outgoing proton mo- 
mentum from the ground-state deuteron wave 
function need not be large, so that the separation 
of proton and neutron can be far greater than in 
the high-Q case. 

The clear implication of this is that for small- 
Q@ stripping at low energies the proton can be well 
removed from the target nucleus when stripping 
occurs, reducing radically thereby the perturbing 
nuclear and Coulomb effects. 

Wilkinson" has pointed out that in such small- 
Q, low-energy reactions, competition through 
the formation of a compound nucleus will in fact 
be discouraged, since the proton which—as has 
been indicated—can be well removed from the 


348 


target nucleus for stripping, will have to pene- 
trate the Coulomb barrier if a compound nucleus 
is to be formed. As in the high-@ case the pro- 
ton must at least partially penetrate the Coulomb 
barrier in order to get sufficiently near to the 
neutron, this effect is much reduced. 

The examples presented in this paper indicate 
an unusually high degree of agreement between 
experimental angular distributions and a simple 
stripping theory uncorrected for Coulomb and 
other effects. A physical description can give a 
clear connection between this agreement and the 
fact that such reactions have low @ values and 
are performed at low bombarding energies. 

A full description and analysis of both experi- 
ments and their results will be published in the 
near future. 

The author wishes to express his gratitude to 
Mr. E. S. Shire for suggesting a study of the 
competition between the stripping and compound- 
nucleus mechanisms for the reaction Li’(d, p)Li®, 
to Professor D. H. Wilkinson for considerable 





guidance and encouragement, and Dr. R. Eden 
for many stimulating discussions. He wouldalso | 
like to place on record his appreciation for finan- 
cial support from the Shell Petroleum Company, 
and from the South African Council for Scientific 
and Industrial Research. 





"This work was done in the Cavendish Laboratory, 
University of Cambridge, Cambridge, England. 

'p. H. Wilkinson (private communication, 1956). 

2S. T. Butler, Phys. Rev. 80, 1095 (1950) 

34. B. Bhatia et al., Phil. Mag. 43 485 (1952); 
Proc. Roy. Soc. (London) A208, 559 (1951). 

‘F. L. Friedman and W. Tobocman, Phys. Rev. 87, 
208 (1952). 

5p. B. Daitch and J. B. French, Phys. Rev. 87, 
900 (1952). 

6 J. Horowitz and A. M. L. Messiah, J. phys. 
radium 66A, 258, 349 (1953). 

™l. P. Grant, Proc. Phys. Soc. (London) A67, 981 
(1954). 

8W. Tobocman, Phys. Rev. 94, 1655 (1954). 

*W. Tobocman and M. H. Kalos, Phys. Rev. 97, 
132 (1955). 

1H. A. Enge and A. Graue, Univ. i Bergen Arbok, 
Naturvitenskap. Rekke No. 13 (1955). 

'D, H. Wilkinson, Phil. Mag. 3, 1185 (1958). 


—_ 





cou 
cou 
Cou 
Sub 


59 


US 


mb 


IS 


— 


81 








VoLUME 3, NUMBER 7 


PHYSICAL REVIEW LETTERS 


OcToBeR 1, 1959 





MEASUREMENT OF u* LIFETIME 


J. Fischer, * B. Leontic, A. Lundby, R. Meunier,t and J. P. Stroot 
CERN, Geneva, Switzerland 
(Received September 10, 1959) 


The conservation of the weak-coupling vector 
current in the presence of strong interactions’ 
is a hypothesis which still lacks convincing posi- 
tive experimental verification. As a consequence 
of this hypothesis a lepton-pion coupling arises 
which is large when the available energy is large 
as in muon capture or short-lived 8 decay. In 
u’ capture the effect amounts to 20-30% of the 
rate, but it seems unfortunately to be quenched 
by an equally large pseudoscalar coupling in- 
duced by virtual 7-meson emission.” Recent ex- 
periments on muon capture between the ground 
states of C’? and B™ support, but do not prove, 
these conclusions.* In 8 decay, relativistic ef- 
fects and radiative corrections tend to eclipse 


the effect of the mentioned lepton- meson coupling. 


The remarkable agreement of the experimental 
with the “theoretical” value for the lifetime of 
the muon (as stressed by Feynman and Gell- 
Mann} is therefore the strongest indication that 
the vector current is conserved and the strengths 
of axial and polar vector currents are equal when 
strong interactions are not present. In view of 
this, we decided, in the spirit of experiments 
involving weak interactions, to attempt to bring 
the experimental value of the yu* lifetime (2.22 
+0.02 psec)* even closer to the “theoretical” 
value (2.31+0.05 psec). Meanwhile, Swanson 
et al.© have reported a measurement where they 
succeeded in this, their new value being 7(:*) 
=2.261+ 0.007 psec. 

The philosophy of our experiment was to ela- 
borate on the detection equipment for muons and 
their decay electrons, and to read their time 
difference in a simple manner from pulses on an 
oscilloscope film. We used 30 hours running 
time of the CERN cyclotron to set up the experi- 
ment and expose the films, while it took us about 
300 hours to read and analyze 10° traces. 

The experimental arrangement is shown in 
Fig. 1. A 180-Mev/c meson beam is momentum- 
analyzed by two successive bending magnets 
after leaving the fringing magnetic field of the 
cyclotron. The velocity-selecting Cerenkov 
counter 1 in coincidence with the scintillation 
counters 2 and 3 registered the passage of a 1. 
Counter 1 consisted of a 10-cm thick cylinder 
subdivided along the beam direction into 5 cells 


which were filled with different mixtures of 
water and ethylene glycol.” The decrease in the 
velocity of the particles by slowing down is thus 
compensated for by an increase in the index of 
refraction (from 1.33 to 1.37) so as to keep the 
angle of the Cerenkov light refracted out of the 
radiator constant. A ring image, corrected for 
spherical, astigmatic, and some of the chromatic 
aberrations, was formed on four 56 AVP photo- 
multipliers by reflecting the light in a nearly 
spherical mirror and passing it through a ring- 
shaped Lucite lens. The detection efficiency of 
the Cerenkov counter for muons was 98 %, while 
for 7 mesons and electrons of the same momen- 
tum it was less than 107°. . - 

The muons were on the average stopped in the 
middle of counter 4, which was a Plexiglas 
cylinder 5 cm thick and 10 cm in diameter. But 
for the side facing the photomultiplier (RCA 
C7071) the radiator had a reflective coating of 
evaporated aluminium, thus making it 47-sensi- 
tive. The counter was surrounded by two layers 
of 1-mm thick mu-metal, which reduced the 
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FIG. 1. Experimental arrangement and a facsimile 


of a typical oscilloscope trace. 
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magnetic field in the interior to of the order of Ntot 
10-? gauss. . 
A coincidence between counters 1, 2, and 3 
with 4 in anticoincidence (1-2-3-4) started the Pe) we oan 


sweep on a 445 Tektronix oscilloscope. The in- 
tensity of the oscilloscope beam was modulated 
with a 10-Mc temperature-controlled crystal 
oscillator. Thus a picture as shown in Fig. 1 
appears on the screen of the oscillograph. The 
pulse from the electron counter 4 in anticoinci- 
dence with 2 (4-2) is shown as a later pulse. The 
accidental counting rate is obtained from the 
frequency of pulses appearing at times earlier 
than the yp pulse, and from the “doubles” rate. 
The latter consists of pulses from a muon, its 
decay electron, and an accidental. We hada 
total of 581 such doubles, and by analyzing their 
time variation along the sweep the corresponding 
variation in the accidental coincidence rate could 
be determined. We considered such an analysis 
necessary since the intensity pattern of the beam 
pulses from a cyclotron can be very asymmetrical, 
but we found a flat distribution over the measured 
time interval. 

Every third trace on the film had an electron 
pulse. The accidental rate was 0.8% of the count- 
ing rate in the first channel. The film was scan- 
ned with the help of a microfilm reader which 
made it possible to locate the position of the pul- 
ses to about 0.02 usec. The scanning was per- 
formed and controlled by five physicists and 
four technicians, and we have not been able to 
discover any systematic errors. The least- 
squares treatment of the data was done with the 
CERN Mercury computer. The results of our 
measurements are shown in Fig. 2. Our value 
for the mean life (2.20+0.015 psec) disagrees 
with the recent result of Swanson et al.,® and 
agrees with the classical value of Bell and 
Hincks.* The error has arbitrarily been in- 
creased from the statistical error (}%). Pre- 
vious experiences with similar measurements of 
lifetimes seem to indicate that errors of less 
than 1% require the convergence of results of 
different experimental methods. 

We gratefully acknowledge the constant experi- 
mental help of P. Duteil and R. Thill and the 
assistance in film reading of P. Meuret and M. 
Namias. 
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GAUGE DEPENDENCE OF THE WAVE-FUNCTION RENORMALIZATION CONSTANT 
IN QUANTUM ELECTRODYNAMICS* 


Kenneth Johnson? 
Physics Department, University of California, Berkeley, California 


and 


Bruno Zumino! 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received August 27, 1959) 


We should like to point out the existence of an 
exact and simple relation between the electron 
Green’s function renormalization constants in the 
general class of “manifestly” covariant gauges. 
We consider those gauges where the unrenormal- 
ized zeroth-order photon Green’s function in 
momentum space has the form 

= 0 
_ ol - bh Ate), (1) 
where A is an arbitrary function of k. By D,,,° 
we shall mean a covariant Green’s function in a 
fixed gauge. 

One can show that the exact, unrenormalized 
electron Green’s functions in the respective 
gauges are connected by the relation 


Gy(x - x’) = expf{ie,?[a(x -x’) - a(0)]}G,(x -x"), (2) 


where e, is the unrenormalized charge of the 
electron, and 


(dk) 


(amyt ee *AlR). (3) 





A(x) = 


This result has been derived previously,’ but 
ina manner which seems to us somewhat ques- 
tionable. An independent and rigorous derivation 
is possible using the techniques of external 
sources, and is given in papers on the gauge 
group in electrodynamics.*** The authors of 
these latter papers were motivated by somewhat 
different considerations. It should also be marked 
that the connection between the Green’s function 
in the radiation gauge [which does not fall into 
class (1)] and the Lorentz gauge has been de- 
tived* by an alternative method which should 
also apply to the case under consideration. 

If the transformation is nonsingular, the func- 
tion A(x) vanishes as x?~«. Further, if the 
singularity of G at p? =-m? is made rigorously a 
pole by cutting off the soft photons with an in- 
variant photon mass yp, then asymptotically in 
any gauge 


G(x -x') = Z,Gm(x -x’) (4) 


as (x -x’)?~«, where 
ip: (x -x’) 


y:p+m 


(dp) e 
(27)* 





Gy (x x!) = ; (5) 
since m, the physical mass of the electron, is 
the lowest invariant frequency contained in 

G(x -x’). Z, is the so-called wave-frequency 
renormalization constant. It follows from Eq. (2) 
that the wave-function renormalization constants 
in the two gauges are connected by the relation 


Z," =exp{ - ie,”d(0)}Z,°. (6) 


To give an example of this relation, one may 
take for D ,, y’ the function 





1 1 
Dy? =p gir aeR): (7) 


where yu is the photon mass and A is inserted so 
that we may regard local electrodynamics as the 
limit of a cutoff theory. Then we choose for A 
the function 


1 1 1 
lorcet @) 
so that we study the one-parameter family of 
gauges for which 


oa 1 1 
Td Be \ (ere) (9) 


In this case, (0) =i()/877)ln(A/y); therefore 








ay =(A/u)"20/?"z,9, 


where a, =e,’/4n is the unrenormalized coupling 
constant. 

If A is allowed to approach infinity, it is clear 
that Z,” can approach a finite and nonvanishing 
limit at most in one of the gauges of class (9). 
Thus, if as A~«, 


(10) 


Z,° ~(A/pL) 


“am, f 
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where zt is finite and not zero, then in the 
gauge in which » =e the renormalization constant 
would be finite and equal to zJ in this limit. 
However, we wish to stress that even if such a 
special gauge did not exist, these considerations 
would in no way affect the possibility that the 
unrenormalized physical theory is consistent. 
Indeed, the reason that a simple relation of the 
type (2) exists is that the gauge change refers to 
the way in which the longitudinal and scalar 
quanta are coupled to the electron, and these 
quanta (because of the conservation of the fluc- 
tuating vacuum currents) enter in an essentially 
noninteracting way. The factor (A/p )¥2o/ 27 is 
produced by the change in the number of these 
quanta in the state in which p* =-m’, and con- 
sequently has nothing to do with the physical 
theory, which is concerned only with the trans- 
verse quanta. Therefore, the fact that these 
degrees of freedom introduce inconsistencies 
only indicates that the use of formulations of 
electrodynamics that employ them is dangerous 
in the investigation of the consistency of the 
(unrenormalized) physical theory at high ener - 
gies. This points to the radiation gauge as use- 
ful for these purposes since it does not fall into 
class (1), unphysical quanta are avoided, and 
the Green’s functions are directly related to 
matrix elements of operators in a Hilbert space 
with a physically acceptable (positive) metric. 

In spite of these difficulties, which have to do 
with the consistency of the local theory, one may 
still make use of the gauges of class (9) ina 
pragmatic way provided that A is chosen larger 
than any energy of physical interest. In this case 
it is useful to note that information about the in- 
frared structure of the Green’s function in var- 
ious gauges may be derived from Eq. (2). It 
should be stressed again that changing the gauge 
alters only the number of scalar and longitudinal 
quanta in the various states and hence has noth- 
ing to do with physics. However, simplifications 
in calculations may occur if the gauge is chosen 
properly. In the low-energy domain (as opposed 
to the infinite-energy domain) the transverse 
quanta should enter in an essentially noninter - 
acting way once the charge is renormalized. It 
is therefore not inconceivable that cancellations 
between real quanta and the longitudinal and 
scalar quanta could be arranged so that a gauge 
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of type (9) might exist where G has rigorously 
a pole at p? =-m’, in the limit as p —0. In this 
case, the singularity in any other gauge can be 
calculated quite simply by using (2). 

If we allow » —0, the difference ie,”[A(x) - d(0)} 
remains finite and becomes equal to 





. H, 2 A(-x?)¥2 ; . = 
ga ft Laldy re) 


Therefore, for A(-x?)”? >1, 
ie,?[A(x) - A(0)] = (ya/27)ln[m(-x?)”?] + const. (11) 


Suppose now that there exists a gauge with 
y =y, where the Green’s function has a pole at 
p? =-m?. Then it is possible to show that in the 
gauge corresponding to an arbitrary value of », 
the Green’s function in the neighborhood of 
p? =-m? will behave like 


(12) 





1 2 Aly -¥)/2m 
yep+m | 
The work of Yennie® indicates that such a gauge 
exists with 
¥ =1-3(a/a,), (13) 


where a is the renormalized coupling constant. 
This is explicitly verified to the fourth order by 
Solov’ev.® 





“This work was supported in part by the National 
Science Foundation and in part by the U. S. Atomic 
Energy Commission. 

tPermanent address: Physics Department, Massa- 
chusetts Institute of Technology, Cambridge, Massa- 
chusetts. 

Permanent address: Physics Department, New York 
University, New York, New York. 
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In this section are printed the abstracts of Articles that 
have been forwarded to The American Institute of Physics 
for publication in THE PHYSICAL REVIEW. In quoting 
information obtained from this section before the appear- 
ance of the corresponding Article, reference should be 
made to “Physical Review (to he published)” rather than 
to this Journal. 


PAIR DISTRIBUTION FUNCTION AND TWO- 
BODY PROPAGATOR. S. Fujita, Department of 
Physics, Northwestern University, Evanston, 
Illinois (Received April 27, 1959). 


Starting from the ordinary definition of the 
pair distribution function, the following theorem 
is proved without the help of the cluster integral 
expansion and “toron” diagrams: The pair dis- 
tribution function in the grand canonical ensemble 
can be expressed in terms of the two-body propa- 
gator corresponding to the scattering in position- 
reciprocal temperature space. A simple expres- 
sion for the internal energy of a nonideal gas 
with pair interactions is obtained in terms of the 
pair distribution function. 


STATISTICAL MECHANICS OF THE STEADY 
STATE. James A. McLennan, Jr., Department 
of Physics, Lehigh University, Bethlehem, Penn- 
sylvania (Received April 28, 1959). 


The description of steady-state phenomena, 
near equilibrium, in terms of Gibbs ensembles 
is discussed. The Liouville equation is modified 
to include external nonconservative forces which 
prevent the system from reaching equilibrium. 
The steady-state ensembles are then obtained as 
(approximately) time-independent solutions to 
the Liouville equation; such ensembles depend 
linearly on the thermodynamic parameters which 
characterize the deviation from equilibrium. 

With the aid of the steady-state ensembles the 
linear relations between the thermodynamic 
fluxes and forces are obtained. 


NUCLEAR SPIN RELAXATION IN LIQUID He’. 
R. H. Romer,* Physics Department, Duke Uni- 
versity, Durham, North Carolina (Received 
April 20, 1959). 


The nuclear spin thermal relaxation time, T,, 
of He* nuclei has been measured between 1.2° and 
3.0° in pure liquid He® under its saturated vapor 


pressure. The measured times are those char- 
acteristic of the bulk liquid. 7, is 300 seconds 
at 1.2° and increases to 550 seconds at 3.0°. 
These results are in fair agreement with the 
theory of Bloembergen, Purcell, and Pound. For 
a 12% solution of He* in He* under its saturated 
vapor pressure, 7, is about 2000 seconds and 
does not show any sharp change at the lambda 
point. In pure He® gas at 4.2°K and 1000 mm 
pressure, T, is at least 1000 seconds. In dilute 
solutions at 1.25°K, 7, is at least 90 minutes for 
a 3.5% solution and at least 120 minutes for a 
1.7% solution. 


. 
On leave from Amherst College, Amherst, Massa- 
chusetts. 


ABUNDANCES OF THE RARE-EARTH NUCLEI 
PRODUCED BY RAPID NEUTRON CAPTURE IN 
SUPERNOVAE. Robert A. Becker* and William A. 
Fowler, Kellogg Radiation Laboratory, California 
Institute of Technology, Pasadena, California 
(Received April 16, 1959). 


Calculations have been carried out, following 
the method of Burbidge, Burbidge, Fowler, and 
Hoyle for the abundances of nuclei in the rare- 
earth region which are produced in the rapid 
neutron-capture process thought to occur in super- 
novae. The recently available rare-earth mass 
differences of Johnson and Bhanot were employed. 
The calculated abundances agree, in general, with 
those given by Suess and Urey. The results of 
the computations support the work of Burbidge 
et al. which showed the effect of spheroidal defor - 
mation above the closed shell at N=126 in enhanc- 
ing the production of Th”*?, U5, U?5*, Cf, etc., 
in supernovae. The effect of different combina- 
tions of temperature and neutron density in en- 
hancing certain relative abundances is discussed 
briefly. 


* 
On leave from the University of Illinois, Urbana, 
Illinois. 


THEORY OF SUPEREXCHANGE. F. Keffer and 
Zs Oguchi, * University of Pittsburgh, Pittsburgh, 
Pennsylvania (Received March 16, 1959). 


The Dirac-Van Vleck-Serber spin operator 
expansion, first applied by Anderson to the 
Kramers superexchange problem, is extended, 
simplified, and systematized in order to handle 
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all overlap contributions arising from a number 
of interacting configurations. The linear cation- 
anion-cation (e.g., Mn**-0~~-Mn**) four - 
electron problem is worked out in detail, taking 
account of all contributions from configurations 
(A) ionic, (B) electron transferred to right, (C) 
electron transferred to left. Group symmetry 
requirements are invoked; and these, together 
with a simple approximation equivalent to 
perturbation theory, are shown to reduce the 
complicated matrix formulation to a single 
linear equation. The solution contains terms 
previously obtained by Anderson, by Anderson 
and Hasegawa, and by Yamashita, and a 
number of important extra terms. All super- 
exchange terms are fourth order or higher in 
the overlap S. A rough numerical evaluation 
with modified Slater wave functions appropriate 
to MnO-type crystals yields an effective super - 
exchange integral of the required size. Brief 
consideration is given to configurations in 
which two electrons are transferred, in parti- 
cular (D) simultaneous transfer of electrons 

to right and to left (Slater mechanism). Unless 
the energy required to form this configuration 
is surprisingly small, its contribution is prob- 
ably not so important, although the problem 
needs to be investigated in detail. Some consid- 
eration is also given to the linear cation-anion- 
anion -cation (e.g. Mn++-Br~-Br~-Mn**) prob- 
lem; the formal solution for the ionic config- 
uration is worked out; and it is shown that super- 
exchange terms first appear in the order S*T°, 
where S is the anion-cation overlap and T is the 
anion-anion overlap. 


“On leave from Tokyo University of Education, 
Tokyo, Japan. 


DENDRITIC GROWTH OF GERMANIUM CRYS- 
TALS. A. I. Bennett and R. L. Longini, West- 
inghouse Research Laboratories, Pittsburgh, 
Pennsylvania (Received January 26, 1959; re- 
vised manuscript received August 13, 1959). 


Controlled dendritic growth of germanium from 
the melt yields long thin strips whose principal 
surfaces are optically flat {111} crystallographic 
planes except for the occasional presence of 
small steps. The crystals grow rapidly in the 
(211) direction, have twin planes parallel to the 
flat surfaces, and are extremely strong mechani- 
cally. The distribution coefficients of impurities 
are close to unity. 
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A mechanism for dendritic growth is proposed, 
in which the presence of at least one properly 
oriented twin plane is fundamental and necessary, 
This mechanism explains most of the observed 
growth features in germanium dendrites, and is 
expected to apply generally to materials with the 
zincblende structure. 

The presence of the twin plane makes growth 


in opposite directions in the twin plane dissimilar, 


not only in the zincblende lattice but generally. 
The effect of this asymmetrization on the growth 
of a-SiC is considered. A crystal growth mech- 
anism based on this asymmetrization is proposed, 
which should be of general validity. 


ENERGY LEVELS FOR RARE EARTH IONS SUB- 
JECT BOTH TO EXCHANGE AND CRYSTALLINE 
FIELDS. Robert L. White and John P. Andelin, 
Jr., Hughes Aircraft Company, Culver City, 
California (Received April 1, 1959). 


The term splitting has been calculated for the 
ground J state of the rare earth ions samarium 
through ytterbium for the case in which the ions 
are simultaneously subjected to exchange and 
crystalline fields of arbitrary relative magnitude. 
The crystalline field is taken to be cubic and as 
representing the potential of eight negative ions 
situated at the corners of a cube. The calculation 
is carried out for the magnetic axis in the two 
major (local) crystalline directions, [100] and 
[111]. The relevance of these calculations to the 
rare earth iron garnets is discussed. 


LOW-TEMPERATURE ULTRASONIC ATTENUA- 
TION IN TIN AND ALUMINUM. Daniel H. Fil- 
son, Department of Physics, University of Cal- 
ifornia, Los Angeles, California (Received 
April 14, 1959). 


According to theory based on an ideal metal, 
the ultrasonic attenuation should be proportional 
to the square of the frequency and to the elec- 
trical conductivity of the metal. Experiments 
were performed to compare theoretical and ex- 
perimental values of attenuation in a frequency 
range from 100 kc/sec to 1 Mc/sec. 

A long wire sample coupled to a barium titanate 
transducer was suspended in a chamber, which 
in turn was placed in a liquid helium bath. A 
short train of sine waves was sent through the 
wire and the amplitude of successive reflections 
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was measured as a function of temperature and 
frequency. The electrical conductivity of the 
sample was also measured. 

Two samples were investigated. High-purity 
tin yielded experimental results in excellent 
agreement with theory. The attenuation in high- 
purity aluminum was proportional to the elec- 
trical conductivity but averaged 45% higher than 
theoretically predicted. 


NEW QUANTUM-MECHANICAL REPRESENTA- 
TION. E. C. McIrvine,* Physics Department, 
Cornell University, Ithaca, New York, and A. W. 
Overhauser,t Ramo -Wooldridge Corporation, 

Los Angeles, California, and Cornell University, 


Ithaca, New York (Received November 17, 1958). 


The superlattice representation is described. 
This representation provides occupation proba- 
bilities which form a well-defined non-negative 
definite analog of the phase space density. The 
superlattice basis functions are orthonormal 
wave packets characterized in general by three 
parameters: a band index, a wave vector, and 
asuperlattice position vector. The crystal 
momentum and position are automatically coarse- 
grained so as to satisfy the Heisenberg uncer - 
tainty principle. The Bloch and Wannier repre- 
sentations are special cases of the superlattice 
representation for particular choices of the 
superlattice parameter. Joint functions of posi- 
tion and momentum, such as the local current 
density, can be represented as the expectation 
value of Hermitian operators through the use of 
this representation. 


*Now at John J ay Hopkins Laboratory, General 
Atomic, San Diego, California. 

‘Now at Scientific Research Laboratory, Ford 
Motor Company, Dearborn, Michigan. 


QUANTUM MECHANICAL APPROACH TO 
THERMAL TRANSPORT PHENOMENA IN 
METALS. E. C. McIrvine,* Physics, Depart- 
ment, Cornell University, Ithaca, New York 
Received November 17, 1958; revised manu- 
script received June 12, 1959). 


A quantum-mechanical transport equation 
‘orresponding to the classical Boltzmann equa- 
ion is developed. Classical transport theory in 
ihe presence of a temperature gradient involves 
the phase-space density, which has no meaning 
Nquantum mechanics due to the complementarity 





of position and momentum. The superlattice 
representation allows the development of a quan- 


tum transport equation from the density-matrix 
Schrédinger equation, through the introduction 

of irreversibility by standard methods. The 
quantum transport equation governing the electron 
occupation probability in the superlattice repre- 
sentation is derived for impurity scattering and 
for phonon scattering. This equation has the 
same field terms as the classical Boltzmann 
equation, but involves the discrete coarse- 
grained wave number and position of the super - 
lattice representation. The scattering terms 
involve transitions between different superlattice 
states, and thus include scattering processes 
which move electrons from one superlattice cell 
to another. The solution of the classical equation 
is affected very little by these quantum effects, 
for all temperature gradients which might rea- 
sonably be encountered. 


* . 
Now at John Jay Hopkins Laboratory, General 
Atomic, San Diego, California. 


DIPOLAR LINE BROADENING AND ENHANCED 
PSEUDO-DIPOLAR MOMENTS. Robert L. White, 
Hughes Research Laboratories, Culver City, 
California (Received May 5, 1959). 


The assumption of a strongly enhanced pseudo- 
dipolar moment in certain ferrimagnetic and anti- 
ferromagnetic materials is shown to be incom- 
patible with the paramagnetic resonance line 
widths in these materials. Why the concept of 
pseudo-dipolar moments is not applicable to line- 
width calculations is discussed. 


ANOMALOUS SKIN EFFECT IN BISMUTH. 
George E. Smith,* Department of Physics and 
Institute for the Study of Metals, University of 
Chicago, Chicago, Illinois (Received April 16, 
1959). 


High-frequency (23.5 kMc/sec) surface resist- 
ance measurements have been made on plane 
surfaces of single-crystal bismuth at 2°K as a 
function of orientation. It has been ascertained 
that extreme anomalous skin effect conditions 
prevail, allowing details of the Fermi surface to 
be deduced from Pippard’s theory. In Shoenberg’s 
model of the electron band, components of the 
inverse effective-mass tensor divided by the 
Fermi energy are found to be a,/E, =9.10, a,/E¢ 
= 0.088, a,/E,=4.7, and a,/E, =0.38 (in units of 
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10°/ev). These results are in essential agree- 
ment with values obtained from de Haas—van 
Alphen experiments and cyclotron resonance. 
The number of ellipses is definitely established 
to be six and the number of electrons found to be 
N=5.5x10'"/cm*. The parameters for the two 
hole ellipsoids are found to be 6,/Ep =8,/Ep, =1.5 
and 6,/E,=0.12. Assuming Shoenberg’s value 
Ee =17.7 Mev, we calculate E, =1.12 Mev from 
specific heat data. It is also found that the re- 
flection of carriers from the surface of the sam- 
ple is predominantly specular in contrast to dif- 
fuse reflection found in other metals. 


. 
Now at Bell Telephone Laboratories, Murray Hill, 
New Jersey. 


STRUCTURE IN THE ENERGY DISTRIBUTION 
OF PHOTOELECTRONS FROM K,Sb AND Cs,Sb. 
E. A. Taft and H. R. Philipp, General Electric 
Research Laboratory, Schenectady, New York 
(Received April 27, 1959). 


The energy distributions of photoelectrons 
from K,Sb and Cs,Sb show structure that is sim- 
ilar in form to structure in the spectral depend- 
ence of the optical absorption. One may ration- 
alize this empirical result by saying that both 
the photoelectric and optical effects arise from 
structure in the state-density of the valence 
band. Assuming that the optical absorption in- 
volves transitions to the conduction band, a 
lower limit for the electron affinity of the crys- 
tal is 0.6 ev for K,Sb and 0.4 ev for Cs,Sb. 


DIFFUSION MEASUREMENTS IN THE SYSTEM 
Cu-Au BY ELASTIC SCATTERING. Robert F. 
Sippel, Field Research Laboratory, Magnolia 
Petroleum Company, Dallas, Texas (Received 
April 24, 1959). 


A broad-range magnetic spectrograph has been 
used to study the diffusion of gold into copper in 
the temperature range 360-500°C by elastic scat- 
tering of protons and deuterons. This is a new 
technique in solid diffusion measurements and 
should have fairly wide applicability in interme- 
tallic diffusion. This investigation extends the 
data a factor of 10° below the sectioning range. 
These results, together with data of sectioning 
observers, show the Arrhenius law to be approx- 
imately valid over a range in D of 5x10°. A 
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precise fit to all the data requires a slight up- 
ward curvature. In the temperature range studig 
here, we obtain Q = 45 750+ 750 cal/mole and D, 
=0.104+0.06 cm?/sec. This represents a de- 
crease of about 4000 cal/mole from the most 
precise data in the sectioning range. This cur- 
vature is attributed to diffusion along internal 
surfaces, although the possibility of multiple- 
diffusion mechanisms still remains. 


ACOUSTIC ATTENUATION IN ALUMINUM DUE 
TO ELECTRON-LATTICE INTERACTION. 
Edward Lax, Department of Physics, University 
of California, Los Angeles, California (Received 
April 22, 1959). 


The attenuation of sound due to the interaction 
of electrons with the lattice has been measured 
in pure polycrystalline aluminum at low tempera- 
tures. Measurements from 26 to 130 kc/sec, 
obtained by recording the decays of the harmonics 
of a longitudinally resonating bar, were com- 
pared with the electrical conductivity of an alu- 
minum wire of identical purity, over a tempera- 
ture range of 3° to 70°K. The attenuation was 
proportional to the square of the frequency. A 
detailed agreement between the shape of the 
attenuation and conductivity curves was obtained, | 
but the experimental attenuation was 50% greater 
than could be explained by the present theory. 


THEORY OF FERROMAGNETIC RESONANCE 
IN RARE EARTH GARNETS. I. g-VALUES. 
C. Kittel, Department of Physics, University of 
California, Berkeley, California, and Hughes 
Research Laboratories, Culver City, California 
(Received April 27, 1959). 


Ferromagnetic resonance phenomena in rare 
earth iron garnets (except gadolinium iron garnet) 
are dominated at room temperature and above 
by the exceedingly rapid relaxation of rare earth 
ions. The rare earth ion relaxation controls the 
g-values and the line widths, as well as the tem- 
perature dependence of these quantities. It is 
shown that in the appropriate limit the g-value of 
the microwave spin resonance line satisfies the 
relation g=g,(M4+Mp)/M,, where g4 and M4 
refer to the ferric lattice and M4, +Mp is the net 
saturation magnetization of the crystal. This 
relation obtains essentially when the B lattice 
relaxation frequency is high in comparison with 
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the AB exchange frequency and relaxation fre- 
quency of the A lattice. The theory accounts 
quite well, with no disposable parameters, for 
the sequence and temperature variation of the 
g-values reported for Dy, Ho, Er, Yb, and Sm 
iron garnets. When at low temperatures the B 
relaxation frequency becomes sufficiently low, 
the g-value should approach the usual result for 
two coupled undamped lattices, in the absence of 
anisotropy effects on the B lattice. The theory 
predicts further that at ordinary temperatures 
the exchange frequency resonance will occur at 
the usual position and its width will be propor- 
tional to the damping constant of the B lattice. 


FERROMAGNETIC SOLUTES IN SUPERCON- 
DUCTORS. Bernd Matthias, Vera B. Compton, 
Harry Suhl, and Ernest Corenzwit, Bell Tele- 
phone Laboratories, Murray Hill, New Jersey 
(Received April 27, 1959). 


Solid solutions of chromium, manganese, iron, 
or cobalt in titanium raise its superconducting 
transition temperature by almost an order of 
magnitude. At very small concentrations this is 
much more than could be expected from a varia- 
tion of electron concentration. 


SATURATION OF THE PARAMAGNETIC RESO- 
NANCE OF A V CENTER. T. G. Castner, Jr.,* 
Physics Department, University of Dlinois, Ur- 
bana, Iliinois (Received April 16, 1959). 


The saturation of the paramagnetic resonance 
of the (halogen), complex (Vi center) in the al- 
kali halides has been studied. The saturation of 
the absorption signal (x’’H,) versus H, is ob- 
tained over a 60-db power range for KCl, KBr, 
and LiF at 78°K. Portis’ theory of inhomogeneous 
saturation has been generalized by omitting the 
assumption that the individual spin packet width 
is very much smaller than the envelope width. 
Methods are developed to determine independ- 
ently, from a given experimental saturation 
curve, the spin packet width 1/T,, the spin-lattice 
relaxation time T,, and the product 7,7,. For 
KCl and LiF values of T,, T,, and T,T, are de- 
termined for the different hyperfine lines of the 
Ve-center spectrum. For KBr, only the product 
T,T, could be obtained. From the results it is 
concluded that the spin packet width is not limited 
by 7,. For KCl at 78°K, T,~77,; for LiF, T, is 





two orders of magnitude or more less than 7, 
and depends in a complicated manner on the ex- 
ternal magnetic field and the angle between the 
V,-center axis and the external field. 

* 


Present address: General Electric Research Labo- 
ratory, Schenectady, New York. 


DEPENDENCE OF THE HOLE IONIZATION 
ENERGY OF IMPERFECTIONS IN CADMIUM 
SULFIDE ON THE IMPURITY CONCENTRATION. 
Richard H. Bube and Arthur B. Dreeben, RCA 
Laboratories, Radio Corporation of America, 
Princeton, New Jersey (Received April 17, 1959). 


The variation of the hole ionization energy of 
imperfections in cadmium sulfide as a function 
of the impurity concentration has been measured 
using photoconductivity in a series of CdS:Ga:Cu 
powders. The Cu concentration varies from 
4x10*" to 2x10” cm™, and each sample was 
prepared with a Cu-concentration to Ga-concen- 
tration ratio of 1.05. The hole ionization energy 
of the sensitizing centers, as determined from 
the thermal quenching of photoconductivity, de- 
creases from about 1.0 ev for low Cu concentra- 
tions to about 0.3 ev for 2x10” Cucm™. The 
results are analogous to other recent findings of 
small hole ionization energies in CdS and CdSe 
crystals. 


MAGNETIZATION OF COMPOUNDS OF RARE 
EARTHS WITH PLATINUM METALS. R. M. 
Bozorth, B. T. Matthias, H. Suhl, E. Corenzwit, 
and D. D. Davis, Bell Telephone Laboratories, 
Murray Hill, New Jersey (Received April 30, 
1959). 


Measurements have been made of the Curie 
points and magnetic moments of three series of 
compounds (Laves phases) MIr,, MOs,, and 
MRu,, where M is a rare-earth element. The 
Curie point is highest (85°K) for the compounds 
containing Gd and falls away when the rare 
earths have larger or smaller atomic numbers. 
This dependence of Curie point on the rare- 
earth element is in accord with the idea that the 
principal exchange interaction is between the 
spins of the 4/ shell of the rare-earth ions and 
the conduction electrons. The magnetic mo- 
ments are closely related to, but differ some- 
what from, the values of Jg for the trivalent 
ions of M. 
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BAND STRUCTURE OF NOBLE METAL ALLOYS: 
OPTICAL ABSORPTION IN a BRASSES AT 4.2°K. 
M. A. Biondi and J. A. Rayne, Westinghouse Re- 
search Laboratories, Pittsburgh, Pennsylvania 
(Received April 22, 1959). 


Calorimetric optical absorption measurements 
at 4.2°K have been made on a representative 
series of a brasses over the wavelength range 
0.23 to 4.0 microns using electropolished bulk 
specimens. Changes in the absorption spectrum 
below 6000 A are interpreted in the light of 
current theories concerning the band structure 
of noble metal alloys. The variation of the infrared 
absorptivity with residual resistivity shows that 
the impurity relaxation time for copper is aniso- 
tropic. This anisotropy increases on alloying and 
suggests that the Fermi surface becomes more 
distorted with increasing solute concentration. 


AGGREGATION OF F CENTERS IN POTASSIUM 
IODIDE. D. H. Goode and P. A. Schroeder, * 
Physics Department, University of Canterbury, 
Christchurch, New Zealand (Received April 23, 
1959). 


As optical aggregation of F centers proceeds 
in KI the R, M, and N bands develop, but later 
lose their identity in a very broad band, which 
is absent in KCl and KBr. The shape and posi- 
tion of this band are temperature independent, 
and it is removed from the peak of the true col- 
loid band. A new band on the short-wavelength 
side of the F band was observed in crystals con- 
taining high concentrations of F centers. 


* 
Now at Division of Applied Chemistry, National Re- 
search Council Laboratories, Ottawa, Canada. 


NATURE OF THE SINGULARITIES IN THE 
SPECTRUM OF A ONE-DIMENSIONAL IONIC 
LATTICE. Joseph Gillis and George Weiss,* 
Department of Applied Mathematics, Weizmann 
Institute of Science, Rehovot, Israel (Received 
April 6, 1959). 


It is known, by a theorem of van Hove, that 
regular lattices will have singularities in their 
frequency spectra, the type of singularity depend- 
ing on the dimension. Although several studies 
of the one-dimensional lattice with Coulomb in- 
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teractions have appeared, no one yet has suc- 
ceeded in identifying the type of singularity that 
appears in addition to the inverse square root 
singularity, in the frequency spectrum. In this 
note we establish that for a one-dimensional 
crystal the singularity is of the form g(w) 
~A{(w-w,)In[1/(w-w,) ]}-*, where A is a constant 
and w, is the position of the singularity. 


“Present address: Institute for Fluid Dynamics and 
Applied Mathematics, University of Maryland, College 
Park, Maryland. 


PROPERTIES OF ALKALI HALIDE CRYSTALS 
AT LOW TEMPERATURES. T. H. K. Barron, 
Department of Physical and Inorganic Chemistry, 
The University, Bristol, England, and J. A. 
Morrison, Division of Pure Chemistry, National 
Research Council, Ottawa, Canada (Received 
March 13, 1959). 


Recent experimental results on low-tempera- 
ture heat capacities and elastic constants of 
alkali halide crystals are examined for evidence 
of anharmonic contributions. It is concluded 
that if such effects are present, their contribu- 
tion to the characteristic temperature at 0°K is 


less than that of uncertainties in the experimental | 


results. 


POSITRON ANNIHILATION IN HELIUM. Thomas 
Bruce Daniel and Robert Stump, Physics De- 
partment, University of Kansas, Lawrence, 
Kansas (Received April 27, 1959). 


Positron mean lives were measured in helium 
gas over a wide temperature range at densities 
from 4.6 to 534 times the STP density and in 
liquid helium at 4.2° and 5.1°K. The orthoposi- 
tronium mean life in the liquid at 5.1°K was found 
to be (9.4+0.6) x10-*® sec which differs little 
from the lifetime in liquid at 4.2°K. Except at 
very low temperatures, the pick-off rate in gas 
agreed moderately well with theory. The low 
pick-off rate associated with orthopositronium 
in liquid was observed also in the gas at temper- 
atures below 9°K for densities greater than 250 
times the STP density. 


"Present address: Midwest Research Institute, 
Kansas City, Missouri. 
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ALPHA-DECAY BARRIER PENETRABILITIES 


WITH AN EXPONENTIAL NUCLEAR POTENTIAL. 


I. ODD-MASS NUCLEI. John O. Rasmussen, 
Lawrence Radiation Laboratory and Department 
of Chemistry, University of California, Berkeley, 
California (Received April 17, 1959). 


Calculation of barrier penetrabilities, reduced 
widths, and hindrance factors for odd-mass 
alpha-particle emitters is made by using the 
diffuse exponential nuclear potential derived from 
optical- model analysis of alpha-particle elastic 
scattering data. The calculations are made on 
the same basis as for even-even alpha emitters, 
reported in a previous publication. 


NUCLEAR BAND STRUCTURE IN Sc“. R. H. 
Davis, Department of Physics, Florida State 
University, Tallahassee, Florida (Received 
April 29, 1959). 


Band structure predictions assuming the forma- 
tion of proton single-particle states above the 
Ca® core in the ground state or one of its ex- 
cited states are compared with the available data 
on the elastic and inelastic scattering of protons 
from Ca®, The band expected above the 3.35- 
Mev state in Ca® is confirmed by experimental 
results, and some evidence is found for bands 
above the higher core states. 


ELASTIC SCATTERING OF 20.35-Mev PROTONS 
BY Zn™, Zn®, AND Zn®*. R. W. Boom” and J. 
Reginald Richardson, Department of Physics, 
University of California, Los Angeles, Califor- 
nia (Received April 27, 1959). 


The absolute differential cross section for the 
elastic scattering of 20.35+0.25 Mev protons has 
been measured for enriched Zn™, Zn®, and Zn®® 
foils. In the angular range of 30°- 160°, about 
50 measurements were made for each foil (spaced 
from 1° to 5°) to an estimated accuracy of about 
5% standard deviation. Scattered protons were 
detected by nuclear emulsions wrapped around a 
4-in. diameter scattering chamber, all angles 
being exposed simultaneously. Detector energy 
resolution is 2.5%, angular resolution is 1° 
standard deviation, and relative angular shifts 
are determined to 0.1°. Correction has been 
made for the finite sizes of beam and detector 
and for multiple scattering in the target and in 
the detector stopper. The +0.25-Mev energy 


spread includes maximum and minimum energies 
due to beam drift, beam spread, and target foil 
thickness. Three minima are found for each iso- 
tope: Zn at 63°, 104°, 142°; Zn® at 62°, 102.5°, 
142°; and Zn® at 61°, 101°, 142°. The absolute 
cross sections are approximately the same ex- 
cept at the third minima, where for Zn™, Zn®, 
and Zn® they are, respectively, 1.46, 1.07, and 
0.61 mb/sterad. 


* 

Now guest physicist at the Institute for Nuclear 
Physics of the Friedrich-Wilhelms University, Bonn, 
Germany. 


DECAY SCHEMES OF THE ISOMERS OF Tc® 
AND Tc*®’. John P. Unik and John O. Rasmussen, 
Lawrence Radiation Laboratory and Department 
of Chemistry, University of California, Berkeley, 
California (Received March 30, 1959). 


An investigation of the decay. schemes of Tc® 
and Tc*” has been made by using high-resolution 
conversion-electron spectrographs, gamma-ray 
scintillation detectors, and coincidence techniques. 

In addition to the 38.9-kev isomeric transition 
in Tc™, eight transitions of the following ener - 
gies have been assigned to Mo*: 204.2, 583.9, 
763, 767.9, 784, 788.0, 822.5, and 837.3 kev. A 
decay scheme is proposed. 

The isomer of Tc®” has been shown to decay by 
a single M4 transition of 96.5+0.1 kev. The ex- 
perimental K:L:M relative conversion-electron 
intensities for this transition are 1:0.48:0.13. 


ANGULAR DISTRIBUTIONS OF PROTONS FROM 
THE (d, p) REACTION WITH DEUTERON ENER- 
GIES BELOW THE COULOMB BARRIER. J. P. 
Schiffer and L. L. Lee, Jr., Argonne National 
Laboratory, Lemont, Illinois (Received April 23, 
1959). 


Cross sections and angular distributions for 
the (d, p) reaction to a known /=1 single-particle 
state of the captured neutron have been measured 
for seven target nuclides between Ti and Ni. 
Targets with thicknesses of several hundred kev 
were used with deuteron bombarding energies 
of 3.8 and 4.5 Mev. The angular distributions 
were found to be similar, with approximately 
2:1 forward peaking and a broad maximum at 
about 60°. In addition, a sharp but relatively 
weak peak was observed at about 25° for the 
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lighter of the target nuclides. This is the angle 
at which the theory of Butler would predict a 
maximum for neutron capture with ‘= 1. Anal- 
ysis of the proton spectra and angular distribu- 
tions indicates that compound nucleus formation 
contributes less than 25% to the reaction yield 
at these deuteron energies. 


BINDING-EFFECT CORRECTIONS IN THE 
ENERGY-LOSS DISTRIBUTION FUNCTIONS 
FOR CHARGED PARTICLES PASSING THROUGH 
THIN ABSORBERS. Walter Rosenzweig, De- 
partment of Radiology, Columbia University, 
New York, New York (Received April 24, 1959). 


In the theoretical treatment of the statistical 
fluctuation in energy loss of charged particles 
passing through thin absorbers, terms which 
are of the order of the electron binding energy 
are ordinarily neglected. The corrections re- 
quired, when this is no longer justified, are 
discussed in particular reference to the develop- 
ment by Symon, which has as its limiting cases 
the Gaussian distribution derived by Bohr and 
Landau’s distribution. Formulas are presented 
which give the correction for each of the three 
weighted parameters appearing in Symon’s deri- 
vation. The corrected parameters can then be 
used in conjunction with Symon’s curves to obtain 
the energy-loss distribution appropriate for a 
given set of conditions. An illustration is pres- 
ented for the cases of 5-Mev protons and 5-Mev 
alpha particles. 


ELASTIC SCATTERING OF N** BY Be®. M. L. 
Halbert and A. Zucker, Oak Ridge National 
Laboratory, Oak Ridge, Tennessee (Received 
April 20, 1959). 


Nitrogen-beryllium elastic scattering was 
measured over an angular range from 32 to 144 
deg in the center-of-mass system with an angu- 
lar resolution of about one degree. The mean 
energy of the incident nitrogen ions was 27.3 
Mev. To distinguish elastic scattering from 
other events both the scattered and the recoil 
particles were detected in coincidence by thin 
CsI(T1) scintillation counters. The elastic scat- 
tering differential cross section is 550 mb/sterad 
at 32 deg c.m. It decreases monotonically and 
more rapidly than csc *(@/2) to a shallow mini- 
mum of about 5 mb/sterad at 106 deg c.m., rises 
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slightly, and then falls to about 2.5 mb/sterad 
at 144 deg c.m., the largest angle measured, 

The data are compared to the predictions of a 
sharp cutoff model for elastic scattering, but 
no agreement is found between this theory and 
the experimental results. 


CROSS SECTION FOR THE Li®°(n, a)H® REAC- 
TION FOR 1.2 <E,, <8.0 Mev. R. B. Murray and 
H. W. Schmitt, Oak Ridge National Laboratory, 
Oak Ridge, Tennessee (Received April 17, 1959). 


The cross section for the reaction Li®(n, a)H* 
has been measured as a function of neutron en- 
ergy in the range 1.2 <E,,<8.0 Mev. An essen- 
tially back-to-back method was used, with a 
thin cylindrical Li®I(Eu) scintillation crystal 
placed concentric with and adjacent to a thin- 
walled ionization chamber containing a deposit 
of fissile material (U*** or Np**’). The magni- 
tude of o,, q 38 measured in this experiment de- 
pends on the absolute value of Ofjgg(U**), while 
the shape of the cross-section vs energy curve 
depends on the known energy dependence of 
Ofjgg(U"** or Np**"). Statistical and other point- 
to-point uncertainties in the data range from +5 
to +9%, while the uncertainty in absolute value 
of the cross section is +7%. The cross section 
obtained from these measurements decreases 
monotonically from a value of 0.28 barn at 
Ey, =1.2 Mev to 0.051 barn at E,, =8.0 Mev. 


ENERGY LEVELS IN A BOUNDED ISOTROPIC 
HARMONIC OSCILLATOR POTENTIAL AND 
NUCLEAR SHELL STRUCTURE. S. Sengupta and 
S. Ghosh, Physics Department, Hooghly Mohsin 
College, Hooghly (W.B.), India (Received March 
14, 1958). 


Energy levels for a three-dimensional bounded 
harmonic oscillator are obtained for intermediate 
distances of the boundary. Energy levels are 
found to depend on the dimensionless parameter 
P,>=(wM/h)R?. pz is the ratio of the strength of 
the oscillator levels (fw) to that of the square 
well levels (t?/MR?). 

Adding a spin-orbit energy 30 times the Thomas 
value, nuclear energy levels are worked out for 
the two mass regions at A =50 and 90. For P=; 
good agreement with the experimental level se- 
quence is obtained. 
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POLARIZATION OF y RAYS IN Pb”. P. H. 
Stelson, W. G. Smith, and F. K. McGowan, Oak 
Ridge National Laboratory, Oak Ridge, Tennes- 
see (Received April 20, 1959). 


Polarization- direction measurements have 
been made on the y-ray transitions in Pb*” re- 
sulting from the decay of Bi?”. These measure- 
ments were made to obtain additional informa- 
tion on both the spin of the level at 2.34 Mev and 
the multipole character of the decay y ray. The 
observed polarization eliminates the possible 
spin assignment 9/2. It is concluded that the 
spin assignment is 7/2 and that the 1.77-Mev y 
ray is predominantly M1. 


TOTAL REACTION AND ELASTIC SCATTERING 
CROSS SECTIONS FOR 22.8-Mev PROTONS ON 
URANIUM ISOTOPES. Clyde B. Fulmer, Oak 
Ridge National Laboratory, Oak Ridge, Tennes- 
see (Received April 16, 1959). 


The total reaction cross sections for 22.8-Mev 
protons on uranium isotopes was determined by 
measuring cross sections for (p, fission), (p,p’), 
(p,d), (p,t), and (p,q@). These data are com- 
bined with previously measured (p,xm) cross 
sections to give total reaction cross sections of 
1.4340.10, 1.44+0.10, and 1.39+0.10 barns for 
22.8-Mev protons on U***, U?*°, and U***, re- 
spectively. Angular distributions of (p,p’), 

(p,d), (b,2), (p,@), and proton elastic scatter- 
ing cross sections are shown for U**> and U*** 


STRANGE -PARTICLE PRODUCTION BY BEVA- 
TRON NEUTRONS IN PROPANE. Charles O. 
Dechand,* Sloane Physics Laboratory, Yale Uni- 
versity, New Haven, Connecticut (Received 

April 24, 1959). 


A liquid propane bubble chamber was exposed 
toa beam of neutrons with energies up to 6 Bev 
from the Bevatron. 10000 pictures of interac- 
tions in the hydrocarbon were scanned to detect 
neutral heavy unstable particles. 349 neutral V 
events were found, most of which came from the 
Stainless steel walls of the chamber. 86% of 
these events could be identified as one or the 
other or either of the neutral strange particles: 
\’ or ©,°. The A°/@,° ratio is about 0.6. 

8200 stars of two or more prongs formed by 
neutrons interacting in the liquid propane were 
observed in the chamber and 17 of these produced 


vs. An additional 5 Vs were formed in single- 
prong events produced by neutrons, and 8 others 
were produced in events in the propane caused 
by charged particles. 

The energy spectrum of the incident neutrons 
was estimated from study of 7-meson production 
interactions in the hydrogen. The distribution _ 
shows that the neutrons had energies up to 6 Bev 
with a mean value of about 4 Bev. For the en- 
ergy range 1 to 6 Bev, the production of strange 
particles occurs in about 1% of all inelastic 
interactions of neutrons with hydrogen and car- 
bon. 


“Now at Combustion Engineering Nuclear Division, 
Windsor, Connecticut. 


K*-DEUTERON SCATTERING IN THE IMPULSE 
APPROXIMATION. Erasmo M. Ferreira,* Im- 
perial College of Science and Technology, Lon- 
don, England (Received September 17, 1959; 
revised manuscript received March 25, 1959). 


It is suggested that a phenomenological analy- 
sis using the impulse approximation of the pro- 
cesses occurring in the scattering of Kt mesons 
by deuterons may be used to get the phase 
shifts for the T=0 isotopic spin state. Typical 
curves are given for the elastic, elastic plus 
inelastic, and charge exchange scattering dif- 
ferential cross sections of 100-Mev K* mesons 
by deuterons on the assumption that only S waves 
contribute and for various ratios of the T=1 and 
T =0 isotopic spin state phase shifts. 


‘On leave of absence from the Centro Brasileiro de 
Pesquisas Fisicas, Rio de Janeiro. 


MEASUREMENT OF PRIMARY DIRECTIONS IN 
EXTENSIVE AIR SHOWERS. C. B. A. McCusker, 
Dublin Institute for Advanced Studies, Dublin, 
Ireland (Received March 24, 1958). 


The abilities of various extensive shower arrays 
to detect possible departures from isotropy of the 
directions of incidence of very high-energy cos- 
mic-ray primaries are analyzed. The effects of 
angular resolution, selectivity of total energy, 
and sensitivity to the nature of the primary are 
discussed. A simple method of analyzing varia- 
tion in intensity with declination is given and 
applied to some recent experimental results. It 
is shown that in two well-known cases these are 


361 








VoLUME 3, NUMBER 7 


PHYSICAL REVIEW LETTERS 











OcToBER 1, 1959 













not consistent with the hypothesis of isotropy of 
incoming directions. 


DIRECTIONAL PROPERTIES OF AN EXTEN- 
SIVE AIR SHOWER ARRAY. C. B. A. McCusker, 
D. E. Page, and R. A. Reid, Dublin Institute for 
Advanced Studies, Dublin, Ireland (Received 
March 24, 1958). 


The directional properties of an extensive air 
shower array consisting of three Geiger-Muller 
counter telescopes have been examined. It is 
shown that the directional properties of this 
array show no great improvement over more 
conventional arrangements. Possible ways of 
improving the performance of the device are 
suggested. 


HEAVY PRIMARY COSMIC RAYS AT GEOMAG- 
NETIC LATITUDE OF 41°N. O. B. Young and 
H. Y. Chen, Southern Illinois University, Car- 
bondale, Illinois (Received April 21, 1959). 


This report includes the results from 9 balloon 
flights at geomagnetic latitude 41°N of altitude 
range from 70000 to 100000 feet. Only prima- 
ries of Z 210 are considered. 2410 tracks are 
involved in Ilford G-5 and G-0 emulsion expo- 
sures. Given are the charge spectra, flux, 
mean free paths, and angular distributions. 


NUCLEAR INTERACTIONS AT ENERGIES OVER 
1000 Bev. J. G. McEwen,* Division of Pure 
Physics, National Research Council, Ottawa, 
Canada (Received April 24, 1959). 


Measurements on two jets, a 0+ 36p jet with 
primary energy about 2500 Bev and a 4+ 29a jet 
with primary energy about 8000 Bev/nucleon, 
are presented. It is shown that the main features 
of the angular distribution of particles from these 
two primary jets and a third energetic secondary 
jet can be explained satisfactorily by the model 
in which mesons created in a nucleon-nucleon 
collision are considered to be radiated isotrop- 
ically from two centers. Further implications 
of this model are discussed; in particular it is 
shown that definite restrictions are imposed on 
the values of energy and angle of emission of 
secondary particles; the average transverse 
momentum is predicted to be relatively constant 
but to have a lower value not only in the forward 
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and backward directions but also around 90° in 
the center-of-mass system of the two colliding 
nucleons. While it is shown that the available 
experimental data are in accord with these pre- 
dictions, more events with energies in the re- 
gion of 1000 Bev must be studied before a definite 
conclusion can be reached. 







































*Present address: Cornell University, Ithaca, New 
New York. 


PRIMARY COSMIC-RAY PROTON AND ALPHA- 
PARTICLE INTENSITIES AND THEIR VARIA- 
TION WITH TIME. Peter Meyer, The Enrico 
Fermi Institute for Nuclear Studies and Depart- 
ment of Physics, The University of Chicago, 
Chicago, Illinois (Received April 22, 1959). 


A series of high-altitude balloon flights was 
carried out in 1957 and 1958 to study the flux of 
primary cosmic-ray protons and a particles 
during variations in the total cosmic ray inten- 
sity. The following results are obtained for a 
particles with energies exceeding 530 Mev/nu- 
cleon under 13.5 g/cm? of air: (a) During a 
large Forbush-type decrease the a-particle and 
proton intensities were closely correlated. This 
demonstrates that a modulation mechanism is 
operating on both components. (b) At certain 
times, variations in the a-particle intensity 
were observed within a few hours which were not 
accompanied by corresponding changes in the 
proton flux. This is tentatively ascribed to an 
anisotropy in the a-particle flux that reaches the 
earth. (c) While there existed an intensity de- 
crease in the proton flux between 1957 and 1958 
which is also observed in the neutron monitor 
station data, no such variation occurred in the 
a-particle flux. A division of the a particles 
into two energy groups (450 Mev/nucleon <E£, 
< 960 Mev/nucleon and E, > 960 Mev/nucleon) 
shows: (a) that the Forbush decrease is of the 
same magnitude in both energy groups, (b) that 
the hourly flux increase observed in some flights 
is about the same in both energy groups, and (c) 
that from 1957 to 1958 the flux in the low-energy 
group increased, while it decreased in the high- 
energy interval, contrary to the well-known be- 
havior of the proton flux. These independent 
a -particle flux variations cannot be explained by 
any of the modulation mechanisms so far proposed. 
We suggest that occasional solar production of 
@ particles may be responsible for our results. 
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The absolute flux of a particles with energies 
exceeding 560 Mev/nucleon at the top of the 
atmosphere was measured on five different days. 


ANALYTIC PROPERTIES OF TRANSITION 
AMPLITUDES IN PERTURBATION THEORY. 
Stanley Mandelstam, Department of Physics, 
University of California, Berkeley, California 
(Received April 17, 1959). 


The analytic properties of two-particle transi- 
tion amplitudes as functions of both energy and 
momentum transfer are examined in perturba- 
tion theory. The modified Nambu representation 
previously proposed by the author for expressing 
these properties is discussed in a little more 
detail. It is shown that, as long as the masses 
do not satisfy certain inequalities connected with 
the existence of anomalous thresholds, the 
fourth-order terms, calculated in the usual 
manner, satisfy the representation. The spec- 
tral functions are calculated explicitly for spin- 
less particles. The proof can be extended to the 
sixth order, but is not worked out here. The 
modifications necessary when there exist anom- 
alous thresholds are mentioned. 


CONSTRUCTION OF THE PERTURBATION 
SERIES FOR TRANSITION AMPLITUDES FROM 
THEIR ANALYTICITY AND UNITARITY PROP- 
ERTIES. Stanley Mandelstam, Department of 
Physics, University of California, Berkeley, 
California (Received April 17, 1959). 


The analyticity properties of transition ampli- 
tudes are used in conjunction with the unitarity 
requirements to generate successive terms in 
the perturbation series, without referring to a 
specific Lagrangian. In the sixth and higher 
orders, production is neglected in the unitarity 
condition; subject to this approximation, it is 
found that the series can be so constructed. No 
analyticity properties which have not been rig- 
orously proved need be employed, and the terms 
are found to satisfy the double dispersion rep- 
resentation. By examining the connection be- 
tween this method and the conventional calcula- 
tion of the perturbation series the types of spec- 
tral function corresponding to different Feynman 
diagrams can be found. Formulas are given for 
the regions in which the spectral functions are 
nonzero. 








TEST OF GLOBAL SYMMETRY IN K-p REAC- 
TIONS. Marc H. Ross and Gordon L. Shaw, 
Indiana University, Bloomington, Indiana (Re- 
ceived April 17, 1959; revised manuscript re- 
ceived July 21, 1959). 


A method of testing the hypothesis that there 
is global symmetry of the pion-baryon inter - 
action is proposed: Upon analyzing low energy 
K~-p scattering data, one finds a variety of scat- 
tering-length solutions which are compatible with 
the elastic scattering, charge-exchange scatter - 
ing, and total charged hyperon production. Our 
suggestion involves the use of the experimental 
=*/Z- ratios to (a) test the global symmetry 
hypothesis (or any other quantitative description 
of the pion-hyperon interaction), (b) reduce the 
ambiguity in the K-N scattering-length solutions, 
(c) hence, predict the £° cross section. We need 
to know the 7 -Y phase shifts, in addition to the 
K-N scattering lengths, in order to predict the 
=*,£~ cross sections. If there is global sym- 
metry of pion-baryon interactions, then we know 
the 7 - Y phases in the absence of a K-N reaction 
channel. We demonstrate how the actual 7 -Y 
phase shifts can be obtained, in a non-perturba- 
tive manner, from the idealized 7 - Y phases 
(i.e., in the absence of K-N reactions) and the 
K-N scattering lengths. Certain fits to the 
present rough data with scattering lengths of 
negative real part are shown to be incompatible 
with global symmetry. Earlier proposals making 
use of the dependence of the hyperon production 
ratios on the 7 -Y phase shifts are also examined 
in terms of this result. The proposed analysis 
involves the assumption that the K is an isotopic 
doublet. 





PROPOSAL FOR DETERMING THE ELECTRO- 
MAGNETIC FORM FACTOR OF THE PION. 
William R. Frazer, Lawrence Radiation Labora- 
tory, University of California, Berkeley, Cali- 
fornia (Received April 23, 1959). 


The possibility of measuring the electromag- 
netic form factor of the pion by extrapolation of 
the cross section for e-+p~-n+n*++e~ has been 
investigated. The method is based on the exist- 
ence of a pole in the pion electroproduction 
scattering amplitude as a function of the invari- 
ant momentum transfer of the nucleon. The 


residue of this pole is the pion form factor mul- 
tiplied by a known coefficient. Since the pole 
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estimate the experimental accuracy necessary 
for a significant extrapolation. Accuracy is re. 
quired which is an order of magnitude better 

than that achieved at present in similar experi- 


ments. 


lies slightly outside the physical region of the 
invariant momentum transfer, an extrapolation 
of the experimental data is required. An ap- 
proximate calculation of the pion electroproduc- 
tion cross section has been made in order to 








